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A Study of the Mechanical Behavior of a 2-D Carbon-Carbon Composite 


(ABSTRACT) 

The objective of this study was to observe and characterize the out-of*plane fracture of 
a 2-D carbon-carbon composite and to gain an understanding of the factors influencing the 
stress distribution in such a laminate. The experimental portion of this study consisted of 
performing an out-of-plane tensile test in a scanning electron microscope and determining the 
modes of failure. Failure was found to be interlaminar, with cracks propagating along the 
fiber-matrix interface. 

Finite element analyses of a two-ply carbon-carbon composite under in-plane, out-of- 
plane, and thermal loading were performed. Stress distributions were studied as a function 
of stacking sequence, undulation aspect ratio, and undulation offset ratio. The results indi- 
cated that under out-of-plane loading o, and t„ were strongly dependent on the geometric 
parameters studied, but o z and 0 , were relatively independent of geometry. Under in-plane 
loading all components of stress were strong functions of the geometry, and large interlaminar 
stresses were predicted in regions of undulation. The thermal analysis predicted the pres- 
ence of large in-.Jane normal stresses throughout the laminate and large interlaminar 
stresses in regions of undulation. 

An elasticity solution was utilized to analyze an orthotropic fiber in an isotropic matrix 
under uniform thermal load. The analysis reveals that the stress distributions in the fiber are 
singular when the radial stillness C r , is greater than the hoop stiffness C M . Conversely, i' 
C„ < C w the maximum stress in the composite is finite and occurs at the fiber-matrr. interface. 
In both cases the stress distributions are radically different than those predicted assuming the 
fiber to be transversely isotropic (C, r = C B0 ). 
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1.0 Introduction and Literature Review 


1.1 Introduction 


Carbon-carbon composites are one of the more unusual composites in that they consist 
of a carbon matrix reinforced with carbon fibers. Since carbon-carbon composites are com- 
posed of elemental carbon, they can withstand very high temperatures (up to 2200 °C) and 
heating rates without appreciable degradation in their properties. In fact, they can actually 
be made stronger at high temperatures. Additional advantages of carbon-carbon composites 
include high thermal shock resistance, high temperature shape stability, chemical inertness 
(except in highly oxidizing environments), and high strength efficiency. Figure 1 compares the 
strength efficiency of carbon-carbon composites with ceramics and superalloys. 1 As 's shown, 
the major advantage of carbon-carbon composites for high-temperature applications is that 
they do not lose strength as the temperature is increased in vacuum. This is in marked con- 
trast to other high temperature materials such as ceramics and superalloys. The figure shows 
two levels of carbon-carbon strength efficiency, .‘he first, labeled Shuttle material, is the 
strength level of the Reinforced Carbon-Carbon (RCC) used on the nose cap and leading 
edges of the Space Shuttle. Although this material is made with low-strength carbon fibers, 


Introduction and Literature Review 


1 



its strength efficiency is still greater than superalloys and ceramics at temperatures greater 
than 1000 °C. More recent research has led to the development of an advanced carbon- 
carbon (ACC), which uses a woven carbon cloth and has twice the strength of RCC and is 
therefore much more strength efficient. 

One of the driving forces behind the development of ACC has been the desire for an al- 
ternative thermal protection system (TPS) for the lower surface of the Shuttle Orbiter. Cur- 
rently the lower surface TPS of the Shuttle consists of High-temperature Reusable Surface 
Insulation (HRSI) tiles. The tiles are a low density silica-based ceramic. Although the HRSI 
tiles adequately withstand the thermal environment (1260 °C) seen by the Orbiter during re- 
entry, the tiles exhibit low strength and lack desirable impact, handling, and other damage 
tolerance characteristics. A thin panel of carbon-carbon, on the order of 0.10 inches thick, 
would be more damage tolerant than HRSi and would most likely result in a significant weight 
saving. Such properties not only make carbon-carbon a desirable material for the TPS of the 
Shuttle Orbiter, but also for more advanced aerospace vehicles. 

Carbon-carbon composites are not without their disadvantages, however. The manufac- 
ture of carbon-carbon composites is a multi-step process, which results in long fabrication and 
processing times and high fabrication costs. The strain-to-faili r e is also low, ranging from 0.3 
to about 1.1%. The high thermal conductivity of carbon-carbon typically requires the use of 
additional insulation material between the carbon-carbon and the aircraft structure. Addi- 
tionally, carbon-carbon composites exhibit low oxidation resistance; therefore, they must be 
treated with oxidation resistant coatings. The general advantages and limitations of carbon- 
carbon composites are have been discussed by Schmidt 2 and are summarized in Table 1. 

It has long been recognized that carbon-carbon composites exhibit low strength in the 
unreinforced directions; therefore, three-dimensionally reinforced carbon-carbon has been the 
norm. 2 3 4 A representative 3-D carbon-carbon is AVCO Mod-3, which has a woven layer of 
orthogonal fibers in the x-y plane and is pierced with fibers in the z direction ’ It has been used 
for re-entry nose tips and for re-entry cassettes, which are containment vessels for radioactive 
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NORMALIZED STRENGTH EFFICIENCY 


STRENGTH EFFICIENCY OF VARIOUS MATERIALS 



Figure 1. Strength efficiency of various materials. 
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isotopes. For rocket nozzle applications cylindrically woven 3-D carbon-car^ on composites 
have been developed and utilized. 

The disadvantages of 3-D carbon-carbon for thin-gauge thermal protection systems are 
four-fold: three-dimensional fabrics are expensive to weave and thus would increase the cost 
of the material; they tend to lack conformability to sharp or multiple contours; parts of difi-orent 
thickness or parts with graded thickness cannot be fabricated; reinforcement in the ou* " >f - 
plane direction reduces the volume fraction of reinforcement in the in-plane directions wit a 
resulting decrease in in-plane properties. 

The logical alternative to 3-D carbon -carbon is 2-D carbon-carbon, which is reinforced in 
the in-plane directions only end thus has none of the disadvantages of the 3-D material listed 
above. However, the absence of thru-thickness reinforcement in 2-D carbon-carbon manifests 
itself in its iow tensile strength in the out-of-plane direction. Current out-of-piane tensile 
strengths are on the order of 50° psi whereas a tensile strength of at le.'ui 1000 psi is needed 
for the foreseeable applications. 

The purpose of this study is to observe and characterize the out-of-plane fracture of 2-D 
carbon-carbon and to gain an understanding of the parameters influencing the out-of-plane 
strength of the material. It is hoped that such an understanding will help guide researchers 
towards a material system that will result in improved out-of-plane tensile properties. 

Chapter 2 presents the results of a;i experimental program in which a carbon-carbon 
laminate was tested under out-of-plane loading conditions. A discussion of the modes of 
laminate failure is given. 

The application of the stress analysis to the region in the carbon-carbon laminate where 
failure occurred is presented in Chapters 3. Chapter 3 focuses on the influence of the woven 
nature of the composite and the ply stacking sequence on the stress distribution in the lami- 
nate. 

Chapter 4 addresses the issue of fiber microstructure, its influence on the transverse 
elastic constants, and the resulting thermal stress distributions. Results are presented for a 
the case of graphite fiber and surrounding matr : x under uniform thermal load. 
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Table 1. General advantages and limitations of carbon-carbon 


ADVANTAGES 

High Temperature Shape Stability 

High Sublimation Temperature 

Low Ablation Recession 

High Strength and Stiffness 
Parallel to Reinforcement Direction 

Retention of Strength and Stiffness 
at Elevated Temperatures 

Thermal Stress Resistant 

Thermal Shock Resistant 

Pseudo-Plastic Mechanical Behavior 

Crack Propagation Resistant 

Nonbrittle Fracture 

Impulse Attenuating 

Chemical Inertness 

Machinability 

Lightweight 

Radi? lion Resistant 

Tailorable Properties 

Nonstrategic Materials 

Fabricability 


LIMITATIONS 

MATERIALS 

Low Off-Axis Mechanical Properties 
Low S’rain-To-Failure 
High Costs 
High Void Content 
Low Fiber/Matrix Bond Strength 
High Thermal Conductivity 
Low Oxidation Resistance 
Low Particle Erosion Resistance 
Nonuniform Pore Distribution 
PROCESSING 

Long Fabrication and Processing Times 
Reproducibility 


DESIGN 

Limited Design and Engineering Properties 
Failure Criteria Lacking 
Attachments and Joints 
Complex Design Methodology 
Complex Response to Environment 
Anisotropic Behavior 

Nondestructive Inspection Not Well Developed 
Limited Applications Experience 
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1.2 Literature Review 


Woven fabric composites have been in commercial use for some time. E'evators on the 
Boeing 727 are a honeycomb structure with woven fabric face sheets. Boeing also uses a 
graphite/epoxy woven fabric composite in the horizontal stabilizer of the 737. McDonnell 
Douglas manufactured DC-10 rudders utilizing woven fabric composites in the early 1970's. 
The wing of the McDonnell Douglas AV-8B Harrier is mostly woven fabric composites. Kevlar 
49 and graphite fabrics are used in the tarings and ailerons, respectively, of the Lockheed 
L-1011. 8 However, research on the mechanical behavior of woven fabric composites has 
generally been neglected by the research community. The exception has been ishikawa and 
Chou, who have published several papers on the in plane elastic behavior of woven fabric 
composites. 7 13 In their analyses, Ishikawa and Chou developed several 2-D models for sim- 
ulating the in-plane behavior of woven fabric composites. In the "mosaic model" 7 8 a plain 
weave composite is modeled as an assemblage of pieces of cross-ply laminates from which 
the upper and lower bounds of elastic constants are estimated. The "fiber undulation 
model" 9 is a refinement of the "mosaic model' that accounts for fiber continuity and the effect 
of fiber undulations on the stiffness of the composite. The mosaic and fiber undulation models, 
however, are insufficient for modeling satin fabrics, where each undulation is surrounded by 
straight fibers. Therefore, a "bridging model" has been developed, which attempts to simulate 
the load transfer between the regions of undulation and the regions of straight fibers. 10 11 
These models have also been extended to the thermal loading case and for nonlinear material 
behavior. 1 * 1 13 

Kriz used the finite element method to study the response of plain weave glass-epoxy 
laminates to cryogenic temperatures. Stiffness and stress distributions were calculated for 
an undamaged composite under thermal-mechanical loading. 14 This was later extended to the 
case of a damaged composite. 15 Free edge effects were also investigated. 15 Kimpara et al. 
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used finite element analysis and acoustic emission to analyze the first knee behavior of a 
fiberglass cloth reinforced plastic. 17 

Jortner presented a mechanistic mode! for the delamination of 2-D carbon-carbon com- 
posites. 18 In this analysis a two-layer model of a plain weave composite was constructed in 
which the weave was approximated by a sine wave. The interface is thus a sine wave, the 
angle of rotation between each point on the interface and the far-field stress was found, and 
the local stresses along the interface were found using standard stress transformation tech- 
niques. Delamination was then predicted using a fracture criteria proposed by Wu and Reuter 
18 that was modified for friction along »he interfaca. 

Stanton and Kipp used the finite element method to model the nonlinear behavior of plain 
weave 2-D carbon-carbon composites. 70 The model was confined to the low modulus material 
used in involute cylinders for rocket nozzle applications. 

Walrath and Adams have done an excellent literature review of the finite element and 
analytical micromechanics and minimechanics models of carbon-carbon composites. 71 They 
review the early efforts to model unidirectional carbon-carbon composites, ncluding attempts 
to examine the influence of the crystallographic orientation of the matrix on the elastic 
moJulus of the composite. Alt-o discussed are the analytical predictions of the stress-strain 
behavior of unidirectional carbon-carbon. The review also reveals that the majority of the 
minimechanics analyses were applied to 3-D carbon-carbon. These analyses ranged from 
relatively simple linear elastic models to those that addressed void content and matrix 
cracking. 
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2.0 Experimental Investigation 


2.1 Introduction 


The experimental phase of this study constituted what could be considered a first step in 
the investiqaiion of the out-of-plane tensile strength of a 2-D carbon-carbon laminate, namely, 
the identification of the failure mode(s). Failure analysis of engineering materials, bn they 
metallic or non-metallic, usually consists of visual and/or microscopic examination of the 
failed components. Investigators don't normally have the luxury of observing failure. If the 
investigators happen to be present, such as in a labora* ->ry environment, they .nay observe 
failure in the macroscopic sense. However, the facilities to observe mic. oscopic events during 
failure aie not usually a /ailable. Thus, the investigator must rely on post-failure analysis and 
engineering judgement in order to determine the cause and mode of failure. The object of this 
experiment was to observe and characterize the out-of-plane tensile failure of a 2-D carbon- 
carbon composite on a microscopic level. This was done by placing a carbon-carbon speci- 
men into a scanning electron microscope (SEM) equipped with a tensile stage and video 
taping an out-of-plane tensile test. This was followed by post-failure microscopy of the frac- 
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ture surfaces of the specimen. The testing was performed at the facilities of the Entry Tech- 
nology Branch, NASA-Ames Research Center, Moffett Field. California. 


2.2 Material Description 


The carbon-carbon panel prepared for this study utilized an eight-harness satin fabric of 
T-300 graphite tows, with each tow insisting of approximately 3000 f:. aments. The fabric 
consists of what is known as warp tows and fill tows. When the fabric is woven the warp tows 
remain straight and the fill tows are arranged such that they go under every eighth tow 
(Figure 2). 

After the fabric is removed from the loom, however, the fabric relaxes a id curvature exists in 
both the warp and fill tows. It should be noted that fabrics do not necessarily have an equal 
number of warp and fill tows. Thus, the mechanical properties in the warp and fill directions 
are not necessarily the same. For example, a T300/5208 graphite cloth/epoxy laminate was 
found to have have an axial re- s of 10 MSI ! n the warp c 1 rection versus 9.4 MSI in the fill 
direct. or. The tensile strengt' ciersed from 83 ksi in the warp direction to 70 k-„ in the fill 
direction. 52 Additionally, fill tows tend to have mc-e twist due to the weaving process, which 
also affects mechanical propert'es. The fabric used in this study had 24 tows per inch in the 
warp direction and 23 tows per inch ir. the f ill direction. The fabric was Fiberite Corporation's 
Style W-133. Its specifications are qiven in Tat 'e 2. 

The panel was supplied by NASA-Langley Research Center. It was prepared by stacking 
nine layers of fabric together in a cross-ply type configuration, with the warp direction of each 
layer oriented 90° with respect to its neighbor. Since there are an odd number of layers, the 
laminate has a warp-predominate direction, a fill-predominate direction, and is therefore un- 
symmetric. The layers of fabric were impregnated with a phenolic resin, cured similarly to 
that of a graphite/epoxy composite, then pyrolized at high temperature in order to convert the 
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Table 2. Graphite Fabric Specifications. 


Fabric Styie 

W-133 

Warp Yam 

T-300 (3K) 

Fill Yam 

T-300 (3K) 

Thickness (mils) 

17.0 1 

(ASTM-D-1777) 

Weight (oz/yd 1 ) 

1080 

(ASTM-D-1910) 

Weave 

8 HARNESS SATIN 

Thread Count; 

Warp x Filling 

24 x 23 

Tracer Yams; 

Warp 

2 INCH CENTERS 

Fill 

6 INCH CENTERS 
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phenolic resin to carbon, it should be not' that the phenolic resin is only about 60 percent 
solids, with the remaining 40 percent volatilizing off. Therefore, it was necessary to expose 
the laminate to several cycles of resin impregnation and pyrolysis in order to obtain a struc- 
turally sound material. The final thickness of the panel was 0.10 inches (2.54 mm). 

Optical photomicrographs of the cross-section of the panel before its initial pyrolysis step 
and after final processing are shown in Figures 3 and 4. respectively. 

Both specimens were prepared for photomicrography by infiltrating them with room- 
temperature hardening resin prior to surface polishing. This step reduces the amount of 
damage done to the specimen during the polishing procedure. Figure 3 shews a partial 
cross-sectional view of the laminate in the as-cured state. This micrograph illustrates the 
woven nature of the composite. Additionally, it shows that there are regions of large vo>ds. 
These voids show up as the dark grey areas as they have been filled in by the resin infiltrated 
into the specimen prior to polishing. They tend to occur at the undulation areas where fill tows 
cross over warp tows. These voids, called macrovoids, occur during the mold step because 
the pressure applied to the laminate is insufficient to fully deform the tows such thai all the 
available volume is filled 

Figure 4 shows a cross-sectional view of the laminate in the fully processed state. The 
lully processed laminate is significantly different from the 3S-cured laminate in that, in addition 
to the fact that the phenolic matrix is now carbonized, the processed laminate contains a large 
number of transverse cracks. Crack density was determined 'rom the micrographs and found 
to be approximately f00 cracks inch. These transverse cracks are most certainly caused by 
the high therma' stresses generated during the high-temperature pyrolysis steps. Addi- 
tionally. there is a significant amount of microporosity, which shows up in the micrograph as 
smail black dots. This is closed porosity, which cannot be removed by the repeated infiltration 
and pyrolysis steps described above 

A low magnification scanning electron micrograph of the composite is shown in 
F'gure 5. This micrograph shows a larger area than Figure 4 and further illustrates the 
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Figiire 3 Optical micrograph of a* cured laminate prior to pyroiyait. 







macroporosity associated with the regions of undulations. Some microporosity is visible, but 
transverse cracks are not. 


2.3 Experimental Procedure 


Several specimens measuring approximately 0.25 x 0.25 in. (6.35 x 6.35 mm) were cut 
from the panel with a diamond saw. The specimens were limited to the above geometry by 
the size and load capacity of the tensile stage. The specimens were then ultrasonically 
cleaned in water in order to remove any foreign matter present due to the cutting process and 
allowed to dry. Each specimen was bonded to two aluminum grips such that its transverse 
cross-section would be the surface of view (Figure 6). The specimens were then sputtered 
with a gold-palladium (Au-Pd) alloy in order to reduce charging effects when in the scanning 
electron microscope (SEM). The specimen to be tested was placed in a Hitachi Model HH-TS2 
tensile stage (Figure 7 and Figure 8a) where the specimen grips were screwed to opposite 
sides of the load frame. The tensile stage was then inserted into a Hitachi Model S-415A 
scanning electron microscope (Figure 8b) equipped with a videotape recording system. 

The videotape recording system was activated and a tensile displacement was applied 
until a crack became visible. The path of the crack was traced until its tip was found. Then 
the load was increased and the propagation of the crack tip along the surface of the specimen 
was recorded. Most often, crack propagation extended beyond the field of view of the micro- 
scope. Thus, it was necessary to go through several steps of applying a displacement, finding 
the new crack tip. increasing the displacement, and video recording crack propagation. Oc- 
casionally. photomicrographs of regions of interest were taken between displacement incre- 
ments. This procedure was carried out until failure of the specimen into two pieces occurred. 
Following failure, one of the two pieces of the failed specimen was remounted in the micro- 
scope such that the fracture surface could be examined and photographed. 
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Figura 5. Scanning electron micrograph of carbon-carbon laminate. 
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Figure 7 . Section view of specimen holder assembly. 
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a) Hitachi Mode! S-415A Scanning Electron Microscope. 
Figure 8. Scanning electron microscope anti tonsils stage. 


b) Hiischi Model HH-TS2 Tensile Stage. 


Reproduced from 
best eveilstxe copy. 
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Fgure 9. Scanning electron micrographs of crack propagation sequence. 
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bers and demonstrates an apparent absence of any fiber-matrix bond. This feature is con- 
sistent with the scahoped crack pattern exhibited in Figure 9d, 

Figure 12, the low magnification micrograph of the fracture surface, shows some regions 
of tow undulations. In the lower-right portion of the micrograph there appears to be a region 
of short fibers (A) aligned 90 0 to the direction of the fibers in the bulk of the fracture surface 
This is not a region of fibers, but a region of matrix similar to that shown at location S in 
Figure 11. This is a result of an undulation from the other half of the failed specimen being 
in contact with it. A similar though less conspicuous region is shown at B in Figure 12. 
However, In this case the undulations are part of the ply shown. 

Intuitively, it would be expected that transverse failure would be interlaminar in nature. 
Intralaminar failure would not be expected because each ply is reinforced n the transverse 
direction due to the woven nature of the fabric. Figure Hand Figure 12 show that the failure 
is indeed interlaminar. The nature of the failure can best be illustrated by considering 
Figure 13. Figure 13 shows a schematic of a two-ply woven fabric laminate. In the particular 
geometric configuration shown, the interface consists o. filaments of similar orientation except 
for the small regions of undulation. A crack propagating along the interface would explain the 
tendency for cracks to propagate between fibers of similar orientation, as exhibited by 
Figure 9a-d, and the morphology of the fracture surfaces exhibited by Figure 11a*b, 

The ultimate out-of-plane tensile strength of the composite was approximately 250 psi 
This differs significantly from the 500 psi tensile strength measured by Maahs ,** It is surmised 
that the difference in strengths is a result of the relative contributions of free edge stresses. 
In the present experiment the specimen was rectangular and much smaller (0.25 in. x 0.25 in.) 
than the round, one-inch diameter specimen used by Maahs, Therefore, the contribution of 
the free edge to the magnitude of the stresses in the interior of Maahs' specimen would be 
expected to to be less than if the smaller specimen were used. Thus, greater strength would 
be exhibited from the larger specimen. 
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Figur# It. Scanning electron micrograph of fracture surface at high magnification. 
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Figure 13. Schematic of location of failure in 2-ply woven fabric composite. 
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2.5 Summary 


It has been shown that the failure of 2-0 carbon-carbon composites due to out-of-plane 
loading is interlaminar. Most significant, however, is the lack of fiber-matrix bond strength, 
as exhibited by the tendency of cracks to propagate along the fiber-matrix interface 
(Figure 10) and the cleanliness of the fiber surface after failure (Figure 11a). Therefore, the 
apparent key to improving the out-of-plane strength of 2-0 carbon-carbon is to improve the 
fiber-matrix bond strength, or reduce the stresses in such regions. 
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3.0 Finite Element Stress Analysis 


3.1 Introduction 


This chapter presents the results of finite element analyses of a two-layer woven fabric 
carbon-carbon laminate under out-of-plane, in-plane, and thermal loading conditions. The 
influence of stacking sequence and fabric geometry on stress states and elastic properties is 
investigated. Of particular interest is the stress distribution at the interface between plies, as 
this is where failure was observed to initiate in experiments on laminates subjected to out- 
of-plane loading (Chapter 5). 


3.1.1 Laminate Geometry 


Consider the geometry of the eight-harness satin fabric exhibited in Fig. 2. Two possible 
stacking arrangements of a two-layer laminate are illustrated in Figure 14 and Figure 15. In 
Figure 14 the laminate is stacked such that the warp tows are adjacent. The cross-sections 
A-A and B-B illustrate the curvatures associated with the warp and fill tows in the laminate. 
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In section A-A the fill tows exhibit curvature and the warp tows do not. In section B-B, how- 
ever, it is the warp tows that exhibit curvature. Borrowing from the nomenclature used to 
describe the stacking sequence of laminates made from unidirectional plies, this stacking 
sequence can be designated as (F/W,/F). For notationai purposes the cross-sections of the 
stacking sequences associated with sections A-A and B-B in Figure 14 are designated 
(F/Wf/F), and ( FIW i IF) w , respectively, where the subscripts f and w identify the tows in the 
cross-section which exhibit curvature. Cross-sections A-A and B-B are identical to those that 
would be exhibited should the laminate be stacked such that the fill tows are adjacent ( i.e 
(W/Fj/W) laminate). A different stacking sequence is illustrated in Figure 15. In this laminate 
the stacking is such that the warp and fill tows alternate. Thus, the stacking sequence may 
be described as (F/W) ? . Comparison of sections A-A and B-B show that one of the sections 
is just the inverted form of the other. 

It should be noted that the undulations are not likely to nest together as shown in 
Figure 14 and Figure 15. A generic offset of the undulations in a (F/Wj/F) w laminate is illus- 
trated i r Figure 16. The offset ratio. R is defined as 

R = X (3.1) 

where 8 and X are defined as shown in Figure 16. Offset ratios of 0, 1. 2, 3, 4, and 5 for a 
( FIW z /F) w laminate are shown schematically in Figure 17. Close examination of the geom- 
etries with R=3 and R = 5 show that one is just the inverted form of the other therefore, of the 
six offset ratios exhibited, only five are unique. Although not pictured, R values of 6, 7. and 8 
would have the same geometries as the inverted forms of R = 2, 1. and 0. respectively. The 
offset ratio R can have values ranging from 0 to 4 in an eight-harness satin. It is important to 
note that R can take on any value between 0 and 4; therefore, there are an infinite number of 
possible offset ratios. It should also be noted that R would be very difficult to control in the 
manufacture of woven fabric laminates. Woven fabrics cannot be handled precisely enough 
to allow laminates with specified values of R to be fabricated using existing fabrication tech- 
niques. 
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Also of interest is the aspect ratio Si of the undulation. £2 is defined as the height H of the 
undulation divided by its wavelength X (Figure 16). The aspect ratio is influenced by both the 
tow size and the processing. A larger tow size, due to either increased filament diameter or 
an increase in the number of filaments in the tow bundle, will result in an increase in Si. In- 
creased autoclave pressure during processing tends to flatten the laminate and decreases 
Si. 

3.1.2 Material Properties 
3.1. 2.1 Matrix Properties 

The matrix material in this study is carbonized Karbon 640, a phenolic resin marketed 
by the Fiberite Corporation. Pyrolysis of the matrix results in a glassy structure, thus the 
matrix can be assumed to be isotropic. The elastic constants of pyrolized Karbon 640 have 
not been published: however, the elastic modulus of pyrolized SC1008, a Monsanto phenolic 
resin with a chemical composition very similar to that of Karbon 640 n , has been determined 
by Bradshaw et al M u to be 4-5 MSI. For this research, an intermediate value of 4.5 MSI was 
assumed for pyrolized Karbon 640. The Poisson's ratio of unreinforced glassy carbon was 
found by Zhao a to lie in the range 0.12-0.17 with the average value being 0.15. Thus, a 
Poisson's ratio of 0.15 is assumed for this study. For an isotropic solid there are only two in- 
dependent elastic constants. The shear modulus G, determined from the above E and v, is 
1.95 MSI. The thermal expansion coefficient of glassy carbon was found by Zhao to range from 
7.9x10 6 /°F to 9.7x10 6 /°F. For this study the intermediate value 8.8x10 6 /°F was used. 
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b) {F/W ? /F) f Stacking Sequence 
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c) (F/W ? /F) w Stacking Sequence 
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Figure 14. (FTWj/F) stacking sequence in a woven fabric laminate. 
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A A 

b) (F/W) 3 Laminate Stacking Sequence A-A 



B B 

c) (F/W) a Laminate Stacking Sequence B-B 


Figure 15. (FA“'1 J stacking sequence in a woven fabric laminate. 
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3.1 ,2.2 Fiber Properties 


The fiber modeled in this study is Thornel T-300. The complete set of elastic properties 
for such a fiber are not available. The difficulty in knowing the fiber properties is complicated 
further by the fact that Eitman 28 estimated that the modulus of a different graphite fiber (T-50) 
increased by as much as 50% when exposed to a 2800 °C carbon-carbon processing cycle. 
Additional work by Maahs 22 on the properties of carbon-carbon suggested that the properties 
of T-300 in a carbon-carbon composite are he same as a virgin T-50 fiber. The axial modulus 
of a T-50 fiber is 57 MSI. 27 Unfortunately, no published data could be found for the other elastic 
properties of the T-50 fiber. Thus, for the purposes of modeling, the remaining elastic prop- 
erties of the T-300 fiber were assumed to be the same as those of Hercules HM fiber because: 
(1) the HM and virgin T-50 fibers exhibit axial moduli which are approximately the same (12% 
difference), and (2) E v v 12 , G„, andv„ values for the HM fiber have been published by Walrath 
and Adams 21 . 

The thermal expansion coefficients were assumed to the same as those of the HM fiber 
and were obtained from Walrath and Adams. They are listed in T jle 3. 

3. 1.2. 3 Tow Properties 

Standard engineering procedure normally entails that the tow properties (in this study a 
'tow' consists of fiber and matrix) be obtained experimentally from a unidirectional carbon- 
carbon laminate. However, the unidirectional material necessary to obtain the elastic con- 
stants is extremely difficult to manufacture. Thus it was necessary to use a micromechanics 
approach and inferences from properties of other material systems available in the literature 
to obtain the tow properties. 

The fiber volume fraction of the tow was determined via optical microscopy to be 0.69. 
Using rule-of-mixtures and the fiber and matrix properties identified in the previous sections, 
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£, of the composite tow was determined to be 40.7 MSI. The remaining elastic constants were 
more difficult to determine due to the unusual nature of carbon-carbon composites. The re- 
sults of the previous chapter indicated a lack of bonding between the fiber and the matrix. 
Therefore, standard micromechanics formulae for the determination of the elastic constants 
of the tow, which assume an ideal bond between the fiber and the matrix, could not be applied. 
Consider, for example, the transverse elastic modulus. The lack of an interfacial bond be- 
tween the fiber and the matrix results in the matrix behaving like a material containing cylin- 
drical voids. The fibers prevent distortion of the holes so the transverse modulus of the tow 
is greater than that of a matrix with unfilled holes. Adams 3 examined the transverse tensile 
and longitudinal shear properties of several unidirectional carbon-carton systems and found 
that the transverse moduli ranged from 0.16-2.1 1 MSI. One system that he studied was the 
T-50/SC1008 system. As previously mentioned the SC1008 resin is virtually the same as the 
Karbon 640 resin used in this study. Adams obtained a transverse modulus of 0.75 MSI for a 
material that had gone through three carbon ization'densification cycles, which is similar to the 
processing that the material in this study went through. Therefore, a value of 0.75 MSI was 
assumed for the transverse modulus E t for this study. 

Adams did not report values of Poisson's ratios or shear moduli, nor has any other ex- 
perimental data for these constants been located in the literature. However, Goetzel w re- 
ported values of Poisson's ratio for a 2-D carbon-carbon material fabricated using a 
thermosetting resin as matrix precursor and a PAN fiber as a fiber precursor. For this com- 
posite v„ z was reported to be in the range 0.10 to 0.17. Since v„ in the 2-D composite is 
somewhat analogous to v 14 in the unidirectional composite, it may be reasonably expected that 
v, 2 will lie in the same range. For this study the intermediate value of 0.15 was assumed. 
Additionally, no data have been located in the literature for v M Recall, however, that the lack 
of a bond between the fiber and the matrix enables the matrix to be modeled as a material 
containing an array of cylindrical voids. Therefore, it can be assumed that the contribution to 
Vj } from the fibers is to prevent distortion of the shape of the voids. Thus, it is expected that 
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v 23 of the tow would vary little from v of the matrix itself. Therefore, a value of 0.15 was as- 


sumed f or . 

If the tow is assumed to be transversely isotropic in the 2-3 plane then the following 
equations hold: 

h 3 = e 2 (3.2) 

v i3 = V 12 (3-3) 


®23 ~ 


2(1 + &v23) 


G 13 - G 


12 


(3.4) 


In order to check if the Poisson ratios used in this study were reasonable, a finite element 
model of a two-layer woven fabric laminate was made in which each layer was approximated 
by a (0/90) laminate. The resulting (0/90), laminate was analyzed and the through-thickness 
Poisson ratio v a was found to be 0.16, which is within the range reported by Gocizel. Thus, 
the values of v used in this study appear to be reasonable. 

Experimentally determined values for the longitudinal shear modulus G, 2 of unidirectional 
carbon-carbon composites are not available. However, since the value of G, 2 of the fiber is 
2.1 MSI, and the value of G for the matrix is 1.9:* MSI, then it can be expected that for an ideally 
bonded composite G l2 would be approximately 2.0 MSI. For this analysis, a value ol 2 0 MSI 
for G, 2 of the tow was assumed. 

The tow thermal expansion coefficients were determined from the fiber and matrix prop- 
erties using the method of Hashin. 29 

The elastic constants of the fiber, matrix, and tow used in this study are summarized in 
Table 3. 
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Table 3. Material Properties. 


Property 

T-300 Fiber 

Glassy Carbon Matrix 

"'vm 

£, (MSI) 

57.0 

4.5 

40.7 

E, (MSI) 

1.5 

4.5 

0.75 

f, (MSI) 

1.5 

4.5 

C.75 

v a 

0.36 

0.15 

0.15 

Vl3 

0.26 

0.15 

01" 

v « 

0.26 

0.1C 

0.15 

G a (MSI) 

0.55 

1.95 

0.326 

G„ (MSI) 

2.1 

1.95 

2.0 

G„ (MSI) 

2.1 

1.95 

2.0 

a, (10 ‘/°F) 

-0.4 

8.8 

-0.074 

a,(10 ‘/°F) 

15.0 

8.8 

12.9 

a,(10 V°F) 

15.0 

8.8 

12.9 
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3.2 Finite Element Model and Implementation 


3.2.1 The Finite Element Model 


As indicated in Figure 14 and Figure 15. a woven fabric laminate is a three dimensional 
structure. In this analysis, however, the laminate is modeled according to the 2-0 cross- 
sectional geometries illustrated in Figure 14b ard c and Figure 15b and c. This is consistent 
with the work o F previous researchers who used two-dimensional analyses to model woven 
fabric lamina in order to avoid the extreme complexity and computational expense associated 
with 3-0 modeling ' * *. It is felt that a 2-0 approach adequately addresses the influence of the 
laminate geometry on the resulting stress state. As such, the warp and fill tows are modeled 
as homogeneous orthotropic sheets, and the laminate is assumed to be homogeneous and 
orthotropic with no v*»r ation in geometry in the y direction. The latter assumption allows a 
generalized plane strain analysis using a two dimensional mesh. Additionally, the and r n 
components of stress are identically zero. 

This model represents the stress state in a generic section A-A or B-fl in Figure 14 and 
Figure 15. As mentioned previously, the geometries associated with sections A-A and B-B in 
Figure 14 are designated the (F.'WjlF), and { F/W r .'F) w laminates, respectively, and the geom- 
etry associated with the cross-sections in Figure 15 is designated the ( F/W) t laminate. Addi- 
tionally. for purposes of analysis and presentation of results, the modeled regions occupy the 
x-z plane. 

One characteristic of a two-layer woven fabric composite is that the geometry is periodic 
in x. Thus it can be divided into repeating units (i.e. 'unit cells'). A typical unit cell of a 
[F/W) t laminate is illustrated in Figure 20a. 

The finite element model is based on the assumptions of linear elastic material behavior 
and linear strain-displacement relations. ups. Eight-node rectangular and six-node triangular 
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plane elasticity, isoparametric elements are incorporated in the model. In addition to these 
six and eight-node elements, a special element has been developed in which the material 
properties are allowed to vary within the element. This eight-node quadrilateral element, 
designated the S8Q. is used to model the curved portions of the tow in the region of undu- 
lation. The finite element model, the concept of generalized plane strain, and the S8Q element 
are fully discussed in Appendix C. 

3.2L2 Finite Element Test Matrix 

The geometric parameters of interest are the stacking sequence, the undulation offset 
ratio R, and the undulation aspect ratio £2. A test matrix was developed which enabled each 
of the above parameters to be investigated. For both in-plane and out-of-plane mechanical 
loading, each of the three stacking sequences. [FIYtJF) w , (FIWJF),. and (F/W) t . were investi- 
gated lor the offset ratios R =0. 0.5. 1. and 4. In each case the aspect ratio £2 = 12 was used. 
This value was obtained by measuring optical micrographs of the material used in the exper- 
imental portion of this study (Chapter 2). 

In order to investigate the influence of £2. additional finite element studies were con- 
ducted for both in-plane ard out-of-plane loading on the (F/Wj/F), laminate for £2 values rang- 
ing from 6 to 24. The offset r jtio R was held constant at 0.5 lor these studies. 

One way woven fabric composites differ from composites made of unidirectional plies is 
that the geometry of woven fabric composites varies with axial coordinate. This is illustrated 
in Figure 18 for two (F/Wj/F). laminates. Each laminate can be divided into two regions of 
primary interest. One of these regions contains the undulation. The other region is desig- 
nated the 'far-field' region, and has the geometry of a cross-ply laminate. Thus, it is expected 
that the e'astic constants and hence the stress distribution in the laminate will vary with axial 
coordinate. In order to investigate the influence of laminate geometry on laminate stiffness, 
average axial stiffness values were calculated for each of the laminates analyzed under in- 
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plane loading. Also, in order to investigate local variations in stiffness within a laminate 
(intralaminate stifTness variations) additional finite element studies were conducted on the 
far-field region and the regions of undulation associated with the (F/kV/F)_ . R~0. 1. and 4 
laminates and the (FIWJF),. R=0 laminate. These regions are illustrated in Figure 19. Ana- 
lyses were performed such that the moduli E„ £,, and £, could be determined. For all cases 
the undulation aspect ratio £2 was held constant at 12. 

As carbon-carbon composites are exposed to elevated temperatures, thermal stress 
distributions in such laminates are very important. Thermal stresses were investigated for a 
[F/WffF), laminate with R=0.5 and £2 = 12. 

The finite element test matrix used in this study is compiled in Table 4. 


3.2.3 Boundary Conditions and Loading 

3.2.3.1 Out-of-Plane Loading 

The boundary and loading conditions for the out-of-plane analysis are illustrated in 
Figure 20b. The lower surface of the laminate (2 = 0) was specified to have zero w displace- 
ment (z direction). A displacement corresponding to an applied out-of-plane strain 
c. = 0.1% was specified at each node on the top surface of the laminate (z = H). The left-hand 
side of the mesh (x=0) was specified to have zero displacement in the x direction. Although 
this boundary condition implies that the mesh has symmetry about the y-z plane, which is true 
only when the offset ratio R equals zero, the results show that this boundary condition does 
not significantly effect the stress distribution in the laminate. The right-hand side of the lami- 
nate (x = L) was free to translate in the x direction, but constrained to remain straight and 
vertical. Absence of this latter boundary condition would have introduced a free edge to the 
model, which does not exist in a unit cell. However, the boundary conditions on the model 
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Table 4. Finite Element Test Matrix. 



Loading j 

In-Plane 

Out-of-Plane 

Thermal 

Laminate 

R 

£2 

R 

Si 

R 

Si 


0.0 

12 

0.0 

H9I 




iiESSii 

12 





(F/Wj/F) w 

1.0 

12 

1.0 

12 

N/A 

N/A 


4.0 

12 

4.0 

12 




0.0 

12 


12 

0.5 

12 



12 

0.5 

12 


- 


1.0 

12 


12 

- 

- 


4.0 

12 

4.0 

12 

- 

- 

(FIYIJF), 

wH 

6 


6 

- 

- 


. 

9 


24 

- 

- 


0.5 

15 

- 

- 

- 

- 


0.5 

18 

* 

• 

- 

• 


0.5 

24 

- 

- 

- 

- 


nm 

12 


12 



(F1W) Z 


12 

in 

12 

N/A 

N/A 


In 

12 

mSm 

12 




4.0 

12 

4.0 

12 
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at x=0 and on the right-hand side (x = L) do require that the interlaminar shear stress i. t 
equals zero at those points. This is discussed in depth later in the text. The final boundary 
condition v,as that the net force in the y direction. F r equal zero. 

1.2.3.2 In-Plane Loading 

The boundary and loading conditions for the in-plane ana'vsis are illustrated in 
Figure 20c. The left-hand side of the mesh (x : =0) was specified to have zero displacement in 
the x direction. In addition, one comer node was specified to have zero displacement in the 
z direction in order to eliminate rigid body displacements. The right hand side was given a 
specified displacement corresponding to an applied axial strain e, = 0.1%. The final boundary 
condition was that the net force in the y direction, F r equal zero. 

3.2.3.3 Thermal Loading 

The boundary and loading conditions for the thermal analysis are illustrated in 
Figure 20d As in the case of in-piane loading, the left side of the mesh (x=0) was specified 
to have zero displacement in the x direction. Also, as in the out-of-plane loading case, the 
right side (x = L) was constrained to remain straight and vertical in order to represent a re- 
peating unit cell. 

3. 2.3. 4 Analysis of Far-Field Region and Regions of Undulation 

The boundary and loading conditions for the analysis of the far-field region and the re- 
gions of undulation are illustrated in Figure 21. The boundary conditions used in order to 
evaluate E, corresponded to those of a uniaxial tensile test (Figure 21a). The left side of the 


Finite Element Stress Analysis 


43 



mesh (x-0) was specified to have zero displacement in the x direction. The right side ( 
x = X) was given a specified displacement corresponding to an applied axial strain 
c, - 0.1%. Additionally, the midplane w displacement at either edge of the mesh (x=0. X ) 
was specified to remain zero. The final boundary condition was that the net force in the y di- 
rection. F y , equal zero. It should be noted that these boundary conditions do not result in a 
state of pure uniaxial tensile stress. Pagano and Halpin (”) s. lowed that significant shear 
stresses are introduced into off-axis tensile specimens due to end grip constraints. The ge- 
ometry of the undulations modeled in this study result in the introduction of z a at the edges 
of the laminate {x=0, k ). 

To evaluate the in-plane stiffness Ey the 'eft side of the mesh (x=0) was constrained to 
have zero displacement in the x direction and the right side of the mesh (x = X) was con- 
strained to remain straight and vertical (Figure 21b). An in-plane force F r corresponding to 
the average in- r r stres o y - 1000 psi was then applied to the model. 

The boundary conditions for determining the out-of-plane modulus E z (Figure 21c) were 
the same as those used for analysis of the full unit cells under out-of-plane loading 
(Figure 20a). 


3.3 Results: Variations in Laminate Stiffness 


3.C 1 Local Variations In Stiffness 


The in-plane elastic modulus E, and the out-of-plane elastic modulus E t of the geometries 
illustrated in Figure 19 were determined by calculating the average stress and dividing by the 
strain. This is expressed mathematically as 
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REGIONS OF UNDULATION 



Figure 18. Regions of interest in (F/W ? /F) w laminates. 
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d) (F/W) v R=0 


Figure 19. Regions of different laminate stiffness in (F/W 2 /F) w laminates. 
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b) O'u-of-Plane Loading, z x . = 0.1% 



c) ln-P!ane Loading, 0 ,. = 0.1% 



Figure 20. Boundary conditions and applied loads for analysis of unit cells it, -* (F/W^/r^ laminate. 


Finite Element Stress Analysis 


47 











(35) 


Ei = 7^A^ idA ' = x,y,z 

The results of this analysis are presented in Table 5. They show that the greatest vari- 
ation in stiffness is in E,. It varies from a maximum of 20.7 MSI in the far-field region 
(Figure 19a) to a minimum of 8.9 MSI in the (F/W% undulation (Figure 19d). This low value 
of E, in the regions of undulation is a direct result of the curvature of the tows. The fill tows, 
which are straight in the far-held region, dominate the axial stiffness. In the undulation, the 
curvature of the fill tows decreases the reinforcement in the x direction and reduces axial 
stiffness (£„). The large difference in E, between the ( F/W) t , R=0 and {F/W^F),,, R=0 undu- 
lations is attributed to the influence of the model geometries or. the deformation of the fill 
tows. In the (F/Wj/F)*, R =0 laminate *he fill tows tend to deform towards the midplane of the 
laminate, but deformation is restricted by the warp tows between the peaks of the undulations. 
The fill tows in the (F/W0 2 laminate, however, are nested such that there is relatively little 
constraint on out-of-plane displacement. Thus, the fill tows in the ( F/W ) z , R=0 laminate are 
more easily able to straighten under tensile load and hence the axial modulus E is reduced. 

The curvature of the fill tows also tends to increase the out-of-plane stiffness E z . However, 
the maximum increase in E t is less than 4 percent. The changes in E x and E z follow the 
transformation equations associated with classical laminated plate theory (CLPT), which state 
that a small variation in fiber orientation from the x axis greatly reduces E , but has negligible 
effect on £,. 

The modulus E y changed very little because E y is dominated by the warp tows, which ex- 
hibit maximum stiffness in the y direction, and whose geometry does not vary with x for the 
generalized plane strain analysis. 
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Table S. Moduli in Specific Regions of 8-Harness Woven Fabric Laminate 


RegiO'-i 

E„ (MSI) 

E y (MSI) 

— 
E t (MSI) 

Far-Field 

20.7 

20.7 

0.77 

(F/Wj/F) wl R=0 

15.2 

2C.4 

0.80 

R = 1,4 

16.3 

20.6 

0.79 

(F/W) 2> R=0 

8.9 

20.4 

0.79 
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3.3.2 Laminate In-Plane Stiffness 


The axial modulus E, was evaluated for the complete repeating unit cell for each of the 
three stacking sequences ((r/l/iyF),, (F/lftyF) w , and {F/W) 2 ) for all four offset ratios (R = 0, 0.5, 
1, and 4). For each of the above cases the undulation aspect ratio fi was held constant at 12. 
Additionally, the influence of Q on E„ was investigated for the (F/Wj/F),, R =0 laminate. Fo r this 
analysis fi values of 6, 9, 12, 18, and 24 were used. 

3.3. 2.1 Effect of Stacking Sequence and Offset Ratio 

The predicted values of E„ as a function of stacking sequence and offset ratio R are pre- 
sented in Table 6. The predicted values of E x range between 16.9 MSI and 19.8 MSI. The 
same maximum E, (19.8 MSI) is predicted for the (F/WJF), and (F/l/vyF) w laminates with R = 0. 
This is attributed to the symmetry' of the laminates. These laminates are symmetric about the 
x-y plane. As discussed in the previous section, symmetry of the laminates reduces the out- 
of-plane deformation that results in lower modulus. Minimum E, is found in the (F!W) V R = 0 
laminate. This result is reasonable because this laminate contains the undulation illustrated 
in Figure 19d, which exhibits the lowest £, of the undulations evaluated (Table 5). 

Comparison of the axial moduli in all three stacking sequences indicate that E, is essen- 
tially independent of offset ratio for R*=0. However, £. appears to be a function of stacking 
sequence. The average values of E, for R ^0 (Table 6) indicate that maximum and minimum 
E, is associated with the (F/W 2 /F)„ and (F/MyF), stacking sequences, respectively. Recalling 
that the (F/VtyF), and (F/W 2 /F) w stacking sequences are associated with cross-sections A-A 
and B-B, respectively, in *he ( F!W Z !F ) laminate (Figure 14), these results indicate that E , and 
E y are not equal in the (F/W 2 /F) laminate, but differ by approximately 10 percent. The average 
value of E k in the (F/W) 2 laminate is 18.4 MSI, which is between those of the 
(F/W^F), and (F/W 2 /F) w stacking sequences. 
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3.22.2 Effect of Undulation Aspect Ratio SI 


The influence of the undulation aspect ratio £2 on the axial modulus E, of the (FIWJF)„ 
R=0.5 laminate is illustrated in Figure 22. The increase in E. frorr. £2=6 to 12=24 is approx- 
imately 6 MSI (44 percent). Most of the increase occurs in the range £2=6-12. where E, in- 
creases from 13.5 MSI to 18.2 MS|. in the range £2=12-24. E, i., creases less than 3%. from 
18.9 MSI to 19.4 MSI. The increase in E, is attributed to the change in curvature of the fill tows. 
The fill tows exhibit less curvature with increasing £2. Thus the axial modulus E, increases 


3.4 Results: Out-of-Plane Loading 


3.4.1 Stress Distributions Ir Regions of Undula. ^n 


As discussed in the previous section, each woven fabric laminate contains a far-field re- 
gion and one or more regions of undulation. The f ar-field reg on behaves essentially like a 
cross-piy laminate, but the response of the undulation region is more complex. To illustrate 
the complex state of stress :n th-. unJulations. stre'.s distributions in ib« {F/WJF), laminates 
are presented in this section 

The stress distributions in the undclations are presented as contour plots. The o. and 
n r components o! stress are discontinuous at the interfaces between the warp and fill tov.*s. 
Additionally, there lends to be large <-.iations in the magnitudes of o, and t„ between the 
warp and fill region: . and are continous at :.ie warp-HI interface). Therefore, the dis- 

tributions m the warp and fill regions are presented in different figures. The contour plots 
presented were generated using tho computer program CAEDS. CAEDS uses the nodal val- 
ues ol the siresses and applies a proprietary algorithm to generate the plot data. It should 
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Table 6. In-Plane Modules E of Woven Fabric Laminate*. 


Laminate 

R 

Q 

E. (MSI) 


0.0 

12 

198 


0.5 


175 

(F/WJF), 


12 

17.5 


4.0 

12 

17.5 

Avg(R*0> 17.5 


wmzam 

12 

19.3 


■ p» 

12 

18.2 

[F/WJF). 


\2 

19.3 


4.0 

12 

19.3 

Avg(R*0) 18.9 


0.0 

12 

16.9 


0.5 

12 

18.5 

{F/W) t 

1.0 

12 

18.3 


4.0 

12 

18.3 

Avg(R*0) 18.4 
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Figure 22. Influence of Undulation Aspect RaUo Li on 


5 ' 



be noted that the algorithm tends to smooth the data, which can result in the loss of areas of 
high stress gradients from the plot. However, within this limitation, the fi ures adequately 
represent the stress distributions in the undulations. 

Contour plots of the normal stresses o, . o r o r and the shear stress for the (F/W/F), 
stacking sequence. R values of 0 and 0.5. are presented in Figure 23 and Figure 25 through 
Figure 31. The contour plots for the cases R = 1 and 4 are presented in Appendix A. 

3 4.1.1 Stress Concentrations 

Prior to presenting the stress distributions in the undulations it is appropriate to address 
the presence of stress concentrations in the laminates. Stress concentrations in structures 
are normally associated with abrupt changes in cross-section. They occur in such frequently 
encountered configurations as holes, notches, and fillets. In heterogenous materials stress 
concentrations also occur in areas where there are point-wise discontinuities in material 
properties. At these point-wise discontinuites the stresses are predicted by linear elasticity 
to be infinite. These points are called 'singularities' or 'singular points', and the stresses are 
said to be 'singular'. Several of the laminates modeled in this study contain singular points 
at the midplane. The number of singular points in each laminate depends on its stacking se- 
quence and offset ratio R. The locations of the singularities m several of the laminates are 
labelled in Figure 18. ,n the finite element method used in this study, the magnitude of a 
stress at a singular point is a function of the density of the mesh at that point. In this analysis 
no attempt was made to refine the mesh near the stress concentration. The density of the 
mesh was constant fcr all cases. Therefore, relative comparisons of the stress distributions 
between different cases can be made. The number and axial coordinate of singular points m 
each of the laminates modeled in this study are listed in Table 7. Their locations are also 
identified as appropriate in the discussions on stress distributions. 
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Table 7. Stress Concentrations in Woven Fabric Laminates 


Laminate 

R 

No. of Singular Points 

Coordinate {-£-) 


0.0 

1 

3.0 

[FWJF), 

0.5 

0 

- 


1.0 

0 



4.0 

0 



0.0 

2 

2.2. 3.8 

(F/WJF). 

0.5 

2 

2.2. 4.8 


1.0 

4 

2.2. 3.8. 4.2. 5.8 


4.0 

4 

2.2, 3.8, 10.2 11.8 


0.0 

0 

- 

(F/W)j 

0.5 

1 

3.0 


1.0 

1 

3.0 


4.0 

1 

3.0 
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3. 4. 1.2 R—0 Laminate 


The contour plot of the in-plane normal stress a, for the {F/WJF )' . R=0 laminate sub- 
jected to out-of-plane loading is shown in Figure 23. Only the upper half of laminate is shown 
since the laminate is symmetric about the x-y plane. 1 he Till tow (Figure 23a) exhibits a fairly 
uniform far-field compressive stress of approximately -100 psi at points A and E. Within the 
undulation, however, there is a significant variation in the stress. At locations where the 
curvature of the fill tow is convex (B. D, and G) the magnitude of o, approaches a tensile stress 
of 3000 psi. Where the curvature of the fill tow is concave (C. F. and H) the magnitude of o, 
approaches -3000 psi. In the warp tow (Figure 23b). the far-field tensile value of o, is 100 psi 
(A and E). It decreases to approximately 70 psi at points F and H and increases to 180 psi at 
point G. The warp tow above the fill tow (B-C-D) exhibits o, values approaching 180 psi at 
points B and 0. but the stress is relatively coi.otant in the remaining area with a value of about 
60 psi. The nonzero far-field o, is a the result of the mismatch in Poisson's ratio between the 
warp and fill tows and the requirement of displacement continuity at the warp-fill interface. 
The warp tow has a higher v a than the fill tow. and thus the fill tows are in compression and 
the warp tows are in tension. The far-field o, in the warp and fill tows are also of equal 
magnitude and opposite sign and thus indicates that equilibrium in both the x and y directions 
is satisfied. 

The stress distributions near the undulation can best be explained by considering the 
deformation of the laminate. Figure 24 shows the deformation of the fill tows in the 
(F/Wj/F), laminates with R values of 0 and 0.5. (The deformations of the other laminates are 
similar ) The fill tows not only undergo displacement in the z direction, but they also contract 
in the x direction and undergo a change in curvature, which is illustrated by tha change in 
angle 0 between the undeformed and deformed fill tows. Thus, the deformation of the fill tow 
is analogous to the bending of a beam. This causes compression along the concave portions 
of the tow (points C. F. and H in Figure 23a) and tension at locations of convex curvature (B. 
0. and G). The axial contraction d in the fill tow (Figure 23a) is indicative of additional 
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compressive 0 , in the undulation. This contraction is responsible for the decrease in a, at 
points F and H in the warp tow (Figure 23b). The increase in a, at point G in the warp tow is 
a result of the bending at that location. It should be noted i!at this laminate does exhibit a 
stress concentation in the form of a singularity at G due to the vanishing warp tows at that 
point. 

The a y stress contour is illustrated in F gure 25. The form of the o } . distribution is similar 
to that of 0 .. but the magnitude of o y is lower. The P'l tow is in tension throughout, with max- 
imum values of o y approaching 140 psi at points B, D, and G and minimum values of 94 psi 
at points C. F. and H. The far-field value of o y is approximately ICO psi. The warp tow is in 
compression throughout, and has a far-field magnitude of approximately -100 psi at points A 
and E. The far-field values of 0 y in the warp and fill tows indicate that equilibrium in the y di- 
rection is satisfied. The normal stress <s y increases to a compressive value of -85 psi at G. 
As with the o, component of stress. a y is a result of the mismatch in Poisson's ratios between 
the warp and fill tows. The small, but definite, variations in a y at the undulation are also a 
result of the deformation of the fill tow. 

The a z stress contour is illustrated in Figure 26. The far-field value of 0 , is approximately 
750 psi in both the warp and fill tows. The maximum values of o, are found in the regions 
F-G and G-H. where the fill tow is inclined and E z is maximum. In these regions c, averages 
800-850 psi with a maximum value of about 1000 psi at points B and D. This distribution is the 
result of a change in the stiffness in the z direction due to the orientation of the fill tow. As 
discussed earlier, the variation in the orientation of the fill tow at the undulation results in 
localized variation in laminate stiffness. In this case the orientation of the fill tow increases 
the stiffness in the z direction Thus, under a uniform applied c, a higher stress is expected. 

The i„ stress contour is shown in Figure 27. The component does not approach the 
far-field value as quickly as the other components of stress. The far-field value of t„ should 
be zero, which corresponds to the value expected from a cross-ply laminate. x„ does reach 
zero, but not within the domain shown in Figure 27. At points A and E, i„ is about 20 psi. The 
largest values of are in the curved regions of the fill tow, between C and D where t, r re- 
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aches a maximum of 300 psi. and between B and C. where x, t reaches a maximum of -300 psi. 
In these locations the orientation of the fill tow with respect to the load axis is analogous to 
that of an off-axis unidirectional laminate under tensile load, and results in a large shear 
stress. 

3.4.1 . 3 R=0.5 Laminates 

The stress contours in the R=0.5 laminate are shown in Figure 28 through Figure 31. 
T he a„ <s y and a, distributions in this laminate are very similar to those of the R=C laminate. 
The shear stresses in the fill tows are also similar. In both laminates the maximum t. in the 
fill tows are ± 300 psi. The major difference between the laminates is in the x n distribution 
in the warp tows. The shear stress in the warp tow region between fill tows (B-C) in the 
R- : 0.5 laminate is approximately 40 psi. On either side (A-B, C-D). i„ is approximately -40 
psi. The difference in between the R=0 and 0.5 laminates is attributed to the offset of the 
undulations in the R=C.5 laminate, which results in greater shear deformation of the layers. 

3.4.1 .4 R—1 and 4 Laminates 

For completeness, the stress contours in the undulations of the R = 1 and R = 4 laminates 
are presented in Figure 99 through Figure 106 in Appendix A. (For the R = 4 laminate the 
stress contour is presented for only one of the regions of undulation. In the R = 4 laminates 
it has been found that the distance between the undulations is large enough that they do not 
influence each o’.her. This is discussed in more detail in the following section on interfacial 
stress distributicns.) Comparison of the figures shows that they all exhibit similar trends in 
their stress distributions. Therefore, the stress distributions are not discussed in detail. 
However, the maximum value of each component of stress for all offset ratios is presented in 
Table 8. 
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A = 500psi 


a) Fill Tow 



A = 1 3 psi 


b) Warp Tows 


Figure 23. Contour plot of 0 < in (F/W 2 /F) ( laminate. R = 0, c* =0.1%. 
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a) Undulation Offset Ratio R = 0 



Undeformed 

Deformed 


b) Undulation Offset Ratio R = 0.5 


Figure 24. Fill tow deformations in (F/W 2 /F) ( laminate. R = 0, 0.5. 1 ^ = 0. 1%. 
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A= 5psi 


a) Fill Tow 



b) Warp Tows 


Figure 25. Contour plot of <T^ in (FW j /F) | laminate. R = 0, =0.1%. 
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The results presented in Table 8 indicates the degree of interaction between the undu- 
lations for all components of stress. The in-plane normal stress o, is greatest for the R=0.5 
case (3300 psi) and decreases with increase of the separation distance between the undu- 
lations (o, = 2690 psi at R=4). The variations in the normal stresses u Y and o, are similar to 
a, in that the maxim _,n values occur in the R =0.5 laminate. The maximum interlaminar shear 
stress in the undulation is relatively independent of offset ratio, differing by less than seven 
percent in the range R=0-4. These results indicate that maximum interaction between the 
undulations occurs for the case R = 0.5. This is attributed to the complex stress state that 
exists because the undulations are close together yet oriented such that the laminate is un- 
symmetric about the x-y plane. 


3.4.2 Interfacial Stress Distributions: (F/W 2 /F), Laminates 


Interfacial stress distributions for the (F/WJF), laminates for R values of 0 and 0.5 are 
presented in Figure 32 through Figure 35. (Here and throughout the remainder of the text, the 
interface is defined as the plane between adjacent woven fabric plies. For the two-layer 
laminate modeled in this study, the interface corresponds to the midplane of the laminate.) 
Stresses were evaluated at the Gauss points nearest the interface for each side of the inter- 
face. In the far-field region the interface is bounded by warp tows on either side. Tne undu- 
lation aspect ratio £} was constant at 1° for all cases. As stated previously, the applied load 
corresponds to out-of-plane strain r. = 0.1% . 

3. 4. 2.1 R—0 Laminate 

In-plane Normal Stresses a > and a : The interfacial a, and c y distributions in the R = 0 iaminate 
are shown in Figure 32 The stress distributions show that the effect of the tow undulations 
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A = 1 10 psi 


a) Fill Tow 



b) Warp Tows 


Figure 27. Contour plot of t in (F/W ? /F) ( laminate. R = 0, =0.1%. 
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A=55Gpsi 


a) Fill Tows 



b) Warp Towr 


Figure 28. Contour ; lot ol o in (FAV *F) ( iariw. ‘e. R=0.5. i. -0.1%. 
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A = i 7 psi 


b) Warp Tows 


Figure 30. Contour plot of O* in (F/W ? /F) ( laminate. R = 0.5, r. 
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b) Warp Tows 


Figure 31. Contour plot of t ^ in (F/W ? /F) ( laminate. R = 0.J 
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•s local. Both components of stress approach far-Hcld values a May from the undulation. The 
distntNition of in- plane normal stress. o„. at the interface is relatively constant in the far-held 
regions at about tOO psi with small decreases at points B and O. which are within the region 
of undulation. A discontinuous jump in o, exists at point C. As discussed in the previous 
section, the presence of a normal stress o, in the far-field regions is due primarily to the 
mismatch in Poisson's ratios between the warp and fill tows. The behavior at points B. C. and 
0 can be explained by recalling the discussion on the deformation of the laminate. The re- 
duction in o. at points B and 0 is due to a compressive component of stress resulting from 
contraction of the fill iow (Figure 24). The jump in o. at point C is a result of the change in 
material properties at that point. Tb® stress at C is evalu^d in a fill tow. which is much stiller 
.i, the x direction than the warp tow. Additionally, the magnitude of o, at C can also be related 
to the deformation of the laminate and the presence of a stress concentration. As discussed 
previously, the geometry at the undulation results in bending of the fill low. This results in 
ar. increase in o. at C. For the plots presented here the stresses are evaluated at the Gauss 
points, which are near but not at t 'oca'.on of ” e stress concentration (C); therefore, the 
contribution of the stress concentration to the magnitude of o. at C is believed to be small 
compared to the contributions from the *.<— in material properties and the bending of the 
tow. 

The interfacial o, d’^fribution (-igure 32b) ■* very similar to that of c, with the exception 
of . iucli lowe.' ma. imum values at C As was the case for o,. the o y stresses are a result of 
the mismatch in Poisson's ratios between the warp and fill tows. The far-field value of o r in 
the warp tows is approximately -110 pM At point C. o r is evaluated in a fill tow, at which the 
magnitude of o t is approximately 130 psi. The decrease in o r at points B and D is a result 
of ihe axial contraction of the undulation (Figure 24) that caused localized decreases at B and 
D in the o. component of stress (Figure 32a) 

Interlaminar Stresses o j and t : The distrbution of the interlaminar normal stress cr, along 
the interface is illustrated in Figure 33a. As was the case for the in-plane normal stresses 
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a, and a r the effect of the undulations is local. Both o t and r a approach far-field values away 
from the undulation. The far-field value of o x is approximately 760 psi (A and E). At points 
B and D, o, decreases to about 740 psi then increases to a maximum of 880 psi at the points 
on either side of C. The decrease in a, at points B and D are attributed to a small compressive 
G t resulting from the change in curvature of the fill tows at the shoulders of the undulation 
(point S). The increase in o, between B and D is due to the increased out-of-plane stiffness ( 
E z ) in the -egion. As discussed previously, the change in orientation of the till tows increases 
the out-of-plane laminate stiffness at the undulation; therefore, under uniform applied strain 
e t , the o t component of stress can be expected to be larger than in the far-field regions. The 
location of the maximum values of o x , on either side of C. appears to be associated with the 
location of largest out-of-plane stiffness along the interface. The local decrease in a z at point 
C can be explained as follows. At point O the curvature of the fill tows is zero; therefore, the 
stacking sequence o? the laminate at C is (F/MyF),.. At this location, the laminate is similar 
to a regular cross-ply laminate. Thus a, has a tendency to approach the far-field value. 

The interlaminar shear stress z„ (Figure 33b.) is essentially 2 ero throughout. This is 
predicted for a cross-ply symmetric laminate. The exception is at C where exhibits values 
of +20 and -20 psi. This is the point of stress concentration, again a result of the idealized 
geometry of the model. 

3 .4.2.2 R=0. 5 Laminate 

In-plane Normal Stresses a and o : The interfacial stress distributions in the R = 0.5 laminate 

» » 

are shown in Figure 34 The distributions of the in-plane normal stresses, a, and o Y are very 
similar to those of the R = 0 laminate. The far-field value of o, (Figure 34a) is approximately 
i r osi and reaches 2100 psi at points C and D. where changes in material properties occur. 
This stress is less than the maximum a, of 2800 psi exhibited by the R =0 laminate This dif- 
ference is due to two factors First, the stiffness of the undulation region changes with offset 
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ratio. Secondly, in the R =0.5 laminate the undulations are ofTset such that no region of stress 
concentration exists. There is virtually no difference in the a y stress distribution (Figu.e 34L) 
between the R=0 and R=0.5 laminates. 

Interlaminar Stresses and y The interlaminar normal st.ess o, (Figure 35a) exhibits a 
far-field magnitude of 760 psi, which is equal to that of the R =0 laminate. Between points B 
and F. o z undergoes changes in magnitude which are caused by changes in stiffness due to 
curvature of the fill tows. The maxirm m value of a z is approximately 840 psi. which is a de- 
crease of about five percent from the R=0 laminate. As with o„ this is attributed to the re- 
duction in stiffness associated with the ofTset of the undulations and the absence of a point of 
stress concentration. The difference of about five percent between the maximum values of 
o z in the R=0 laminate which has a singular point, and the R-0.5 laminate, which lacks a 
singular point, suggests that for this analysis the influence of the singularity on the stress 
distribution in these laminates is minimal. As discussed previously, this may be a result of 
the density of the finite element mesh at the singular points. Had the mesh been finer at the 
singular points their effect might be more noticeable. 

The distribution of in the R = 0 5 laminate (Figure 35b) differs significantly from that of 
the R = 0 laminate. Whereas the R=0 laminate exhibited essentially zero shear stresses, t iz 
in the R=0.5 laminate ranges from -50 psi at points B and F to 40 psi at point D. This is a 
result of the offset of the unaulations in the laminate. As discus sed earlier, the fill tows un- 
dergo a decrease in radius of curvature during deformation. This results in greater con- 
traction of the fill tows in the x direction than is attributed to standard Poisson ratio effects. 
In this laminate, the offset of the undulations causes a mismatch in the contraction between 
the two plies of the laminate. Therefore, a significant shear stress is generated at the inter- 
face. 
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3. 4. 2. 3 R—1 end 4 Laminates 


The interfacial stress distributions in the R = 1 and R = 4 laminates are shown in 
Figure 115 through Figure 118 in Appendix A. The distributions of the normal stresses o„ 
a y , a, in the laminates differ very little from the R = 0 and R=0.5 laminates, and their further 
discussion is not warranted. (The differences and similarities between the laminates are 
summarized in a later section.) However, significant differences do exist for the interlaminar 
shear stress t„. The distribution in the R = 1 laminate is shown in Figure 116b. In this 
laminate, t„ ranges from approximately 12 psi to -48 psi. In the R=4 laminate (Figure 118b) 
t„ ranges from 30 to -30 psi. Thus, the magnitude of decreases with separation distance 
(R) between the undulations. This indicates that there is significant interaction between the 
undulations, with respect to the component of stress, for small non-zero offset ratios. 

Of special interest is the shear stress distribution in the regions -j- = 0-2 and 
-j- = 12-14 in the R=4 laminate (Figure 118b). In the large far-field region between points B 
and C, is non-zero, ranging from -7 to -15 psi. This implies that should also be non-zero 
at A and D This contradicts the boundary condition of the finite element model which requires 
x„ = 0 at the boundaries of the unit cell. Clearly then, modeling a single unit cell does not 
accurately predict the the distribution in the R = 4 laminates. To address the question of 
the accuracy of the one unit cell model, a two unit cell model was generated and the results 
from the two unit cells were compared. This v‘ as done for ail components of stress for the 
(F/W 2 /F), , R = 4 laminate for both m-plane and out of-plane loading. The results, which are 
presented in detail in Appendix 3, indicate that the only deviation in the results between the 
one and two unit cell models is in the x, t for the out-of-plane loading condition. Still, there 
was no difference between the two models in the maximum predicted vclje of t„. Thus, the 
one unit cell model is considered sufficient for modeling the t„ distribution for all offset ratios. 
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3.4.3 Interfacial Stress Distributions: (F/W y /F) w Laminates 


The interfacial stress distributions for the (F/W t fF) w laminates are shown in Figure 36 
through Figure 43. These laminates are different from the \FI\NJF), laminates in that the warp 
tows exhibit curvature (Figure 14c). Additionally, the geometry of the (FIWjlF) w laminate is 
such that it contains more sources of stress concentration than the {F/WJF), laminates 
(Table 7). 

3. 4. 3.1 R-0 Laminate 

ln*plane Normal Stresses o x and a : The a, and a y stress distributions in the R = 0 laminate 
are shown in Figure 36. The far-field o, and o y values are -110 psi and 110 psi. respectively. 
Both components of stress are maximum at points B and D. In this case, the high values of 
stress at these points are due to bending of the warp tows. The decrease in a, and o y at C is 
a result of the change in material properties at the interface between B and D. In this region 
both components of stress are are evaluated along a fill tow, which is in tension in the x di- 
rection and compression in the y direction due to Poisson effects. 

Interlaminar Stresses o ; and x : The distrl jtions of a z and x„ are shown in Figure 37. The 
magnitude of c, is 760 psi in the far-field regions while the maximum values of approximately 
950 psi occur at points B' and O'. These points correspond to the locations along the interface 
of the laminate where the maximum values of E t exist. At C the geometry of the laminate is 
similar to that of a cross-ply laminate. Thus a, tends toward the far-field value. In the far-field 
region is zero, which is consistant with the response of a cross-ply laminate. This laminate 
is symmetric, so x„ would be expected to vanish along the interface. The magnitude of at 
points B' and D' and the presence oi oscillations in the stresses near these points, are due to 
numerical instabilities associated with the presence of stress concentrations in the model. 
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Figure 35. Interfacial o ^ and t 
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3.4.32 R—0.5 Laminate 


In-plane Normal Stresses o ( and The 0 . and O r stress distributions in the R=0.5 laminate 
are shown in Figure 38. The high values Of a t at points B. D. E. and G. relative to the far-held 
o„ are due to bending of the warp tows. The presence of negative components of o. at points 
C and F requires some discussion. Consider the negative a, stress at point C. which occurs 
in the warp tow below the midpiane of the laminate. The presence of this compressive stress 
is caused by the deformation of the curved portion of the warp tow in the top layer. Defor- 
mation of the warp tow causes its shoulders (points C and E) to be pulled in toward each other. 
The result is compression in the warp tow below the interface at point C. The response at 
point F is similar however, in this case the deformation of the bottom layer fill tow causes 
compression in the top layer fill tow at F. The a, distribution is very similar to that of the 
R=0 laminate. The peaks in a r at points B. D. E, and G are due to bending of the fill twos at 
the undulations. The negative components of a y are due to changes in material properties 
along the interface. 

Interlaminar Stresses a and t : The interfacial stresses o. and x„ for the R =0.5 laminate are 
shown in Figure 39. The a t is very similar to that of the R=Q laminate. The interlaminar 
shear stress x a is non-zero due to the asymmetry of this laminate. Points B and G are sources 
of stress concentration. 

3. 4. 3. 3 R—1 and 4 Laminates 

The interfacial stress distributions for the R = 1 and R=4 laminates are shown in 
Figure 40 through Figure 43. These laminates contain four sources of stress concentration, 
compared to only two each for the R = 0, 0.5 laminates. This is especially noticeable in the 
o, and t„ distributions (Figure 41b. and Figure 43b). where continuity of these components 
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of stress does not exist at the sources of stress concentre ~ Generally, the stress distrib- 
utions in the R = 1. 4 laminates are very similar to those of t. -0. 0.5 laminates. The sim- 

ilarities and differences between them are summarized in section 3.7. 


3.4.4 Interfacia! Stress Distributions: (F/W), Laminates 


The stacking sequence of this laminate differs from the other laminates considered in that 
it is not symmetric in the far-field region. This results in the midplane being bounded by tows 
of opposite orientations. The tow above the interface is a warp tow and the tow below the 
interface is a fill tow. 

3. 4.4.1 R=0 Laminate 

In-plane Normal Stresses o < and a : The interfacial o x and o y distributions for the R=0 lami- 
nate are shown in Figure 44. The distributions in this laminate are very similar to the dis- 
tributions in the laminates previously discussed. This laminate exhibits far-field o„ and o y 
values of 110 psi and -110 psi, respectively, in the wa-p tow (above the interface), with a dis- 
continuity at C due to a change in material properties. The fill tow (below the interface) has 
far-field a, and a y values of -110 psi and -110 psi, respectively, with large changes in magni- 
tude at points B and 0 The change in a, at B and D is due to bending of ' le fill tows at those 
points. The change in o y at B and 0 is due to a change in material properties. Between B 
and D, o y is evaluated in a warp tow. 

Interlaminar Stresses a z and t m : The o, and x„ distributions (Figure 45) are very simifar to 
those of the other laminates. Maximum values of a, occur at B and D where the out-of-plane 
stiffening effect of the undulations is greatest. The interlaminar shear stress x w shows maxi- 
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Figure 42. Interfacial and in (F/Wj/F)^ laminate. R=4, e j *0.1%. 
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mum values at points B and 0 (Figure 45b). It is at these points where maximum interfacial 
shear exists. Of particular interest is the difference between the two x„ distributions pre- 
sented in Figure 45b. Satisfaction of the stress continuity requirements at the interface re- 
quires that t m be the same in both layers. In this case, continuity of x„ is not satisfied between 
points A and B and between D and E. This is due to the presence of a coarse finite element 
mesh in these regions. The mesh density between B and D was four times greater than be- 
tween A and B and D and E. The existing mesh was fine enough to satisfy continuity at the 
points of maximum shear stress (B and D). 

3. 4. 4. 2 R=0.5 , 1 and 4 Laminates 

The remaining laminates in this stacking sequence (those with R=0.5, 1 and 4) exhibit 
stress distributions similar to the previously discussed; therefore, the specifics of the distrib- 
ution of each component of stress for each of the remaining laminates will not be presented. 
For completeness, however, the interfacial stress distributions are presented ! n Appendix A. 
A summary of the results for out-of-plane loading is presented in section 3.4. 


3.4.5 Influence of Undulation Aspect Ratio £2 


The influence of the undulation aspect ratio (£2) on the maximum interfacial stress was 
evaluated for a (h/W 2 /F), laminate with an offset ratio R=0.5. The maximum values of o„ o y , 
Oj, and along the interface, as a function of £2 are shown in Figure 46 through Figure 49. 
The in-plane normal stress a, exhibits the greatest dependence on £2 (Figure 46). At £2 = 6, 
o, is approximately 3200 psi. It decreases to 2100 psi at £2 = 12 and to approximately 700 psi 
at £2 = 24. The decrease in maximum c, is due to the reduction in curvature of the fill tow with 
increasing £2. As di scussed previously, the deformation of the fill tow is analogous to a beam 
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Figure 44. Interfacial o and a In (F/W). laminate. R * 
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bending problem. The increase in £1 effectively 'straightens' the fill tow and thus reduces the 
a M component of stress attributed to bending. 

The other components of stress, a y . a „ ana are effected less by £2. The decrease in 
o y (Figure 47) is less than 10% between £2=6 and $2=12, and does not decrease for 
$2 > 12. The interlaminar normal stress a z (Figure 48) decreases from approximately 1000 
psi at £2=6 to 820 psi at £2=12, and is approximately 780 psi at £2=24. The interlaminar shear 
stress t„ (Figure 49) decreases from 78 psi at £2 = 6 to 52 psi at £2=12 and equals psi 40 at 
£2=24. 
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Ficiire 46. Maximum interfacial O as a function of £}. r. =0.1%. 
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Figure 48. Maximum interfacial O f as a function of r.^ — 0.t%. 
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Figure 49. Maximum interfacial t as a function of Q. =0.1%. 


Finite Element Stress Analysis 


97 



3.6 Results: In-Plane Loading 


3.5.1 Laminate Deformations 


As was the case for out-of-plane loading, the response of the woven fabric laminates to 
in-plane loading (e, = 0.1%) is best understood by first considering the deformations of the 
laminates. The deformations of the (F/WJF), laminates for R values of 0 and 0.5 are illustrated 
in Figure 50. The R=0 laminate (Figure 50a-b) illustrates the tendency of the curved tow to 
straighten under axial loading. As was the case for out-of-plane loading, the deformation of 
the fill tows is analogous to the bending of a curved beam. Figure 50c-d illustrates the effect 
of a lack of symmetry on the deformation of the R=0.5 laminate. The laminate is straight in 
the far-field regions (which is expected because the laminate is symmetric in those regions) 
and curved in the region of undulation. The deformation is similar to that of a single lap joint. 

The deformations of the fill tows themselves are illustrated in more detail in Figure 51. 
The fill tow in the R=0 laminate (Figure 51a) shows not only translation and axial extension, 
but also a reduction in curvature. This is in contrast to the out-of-plane loading condition, 
where the tendency is for the curvature of the fill tows to increase. The hi; tows of the R =0.5 
laminate (Figure 51b) exhibits the same characteristics of the R=0 laminate. Additionally, 
there is cantilever-type bending of the tows. As discussed previously, this is due to the lack 
of symmetry of the laminate. 
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a) Undulation Offset Ratio R=0, undeformed 



b) Undulation Offset Ratio R=0, deformed 



c) Undulation Offset Ratio R =0.5, undeformed 



d) Undulation Offset Ratio R =0.5, deformed 


CD 


Figure 50. Deformations of (F/W 2 /F) f laminates. 6^= 0.1%. R = 0, 0.5. 





Undeformed 

Deformed 


a) Undulation Offset Ratio R= 0 
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Deformed 


b) Undulation Offset Ratio R=0.5 


Figure 51. Fill tow deformations in (F/W 2 /F) f laminate, 0.1%. R =0,0.5. 
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3.5.2 Stress Distributions in Regions of Undulation 


In this section the stress distributions in the undulations are discussed for the in-plane 
loading condition. Results are presented for the [F/Wf/F), laminates for R values of 0, 0.5. 1, 
and 4. 


3. 5. 2.1 R=0 Laminate 

The in-plane normal stress o x for the R=0 laminate subjected to in-plane loading is 
shown in Figure 52. This laminate is symmetric about the x-y plane, thus only the results for 
the upper half of the laminate are presented. The laminate is in axial tension throughout, with 
far-field o x values (locations A and E) of approximately 40 ksi in the fill tow and and less than 
1 ksi in the warp tows. The fill tow exhibits regions of local maximums at locations C, F and 
G, where o x is approximately 54 ksi. The minimum values of o x in the fill tow are at locations 
B and D, where a, it is approximately 20 ksi. A local minimum of 34 ksi exists at location G. 
The overall stress distribution is a result of contributions from the bending of the tow as well 
as the axial extension of the tow. In this case, the contribution of the bending of the tow is to 
reduce c, in the fill tow at locations B, D, and G, and to increase o, at locations C. F, and H. 
In the warp tow, maximum deformation occurs at point I, and is caused by the deformation 
of the fill tow below it (Figure 50a). The deformation of the fill tow causes the warp tow above 
it to bend as it is translates in the z direction, thus increasing o, to approximately 1.5 ksi at I. 
A local maximum in o, also exists at location G. The reason for this is two-fold. First, the 
reduction in axial stiffness in the region of undulation, which is duo to the curvature of the fill 
tow, results in a larger component of the load being absorbed by the warp tows. Secondly, 
location G is a point of stress concentration in this laminate (Table 7). 

The a y stress contour is illustrated in Figure 53. As was the case for out-of-plane loading, 
the presence of a component of stress is due to the mismatch in Poisson's ratios between 
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the warp and fill tows. In the Till tow c y is maximum at location G, reaching a value of ap- 
proximately 850 psi. Local maxima of approximately 200 psi are found at locations F and H. 
The far-field value of a y (locations A and E) is about 100 psi. In the warp region the maximum 
value of o y is 740 psi at location G. A local maximum of 50 psi is found at location I. The re- 
mainder of the warp region is in compression, and is about -100 psi. The far-field values of 
o y in the warp and fill tows indicate that equilibrium in the y directions is satisfied. 

The o z stress contour is shown in Figure 54. The interlaminar stress o, is zero in the 
far-field regions, which is consistent with CLPT. In the undulation a z is nonzero and reaches 
a maximum value of approximately 5 ksi at location G. The magnitude of o z resulting from 
in-plane loading is significantly greater than that due to out-of-plane loading. The maximum 
value of a t under out-of-plane loading for this laminate was less than 1 ksi (Figure 26). 

The stress contour is shown in Figure 55. The interlaminar shear stress x„ is zero in 
the far-field regions (locations A and E). but varies from -11 ksi to 11 ksi in the curved portion 
of the laminate between locations F and H. As was the case for out-of-plane loading, the ori- 
entation of the fill tow results ii. a response that is analogous to that of an off-axis 
unidirectional coupon under tensile load. The orientation of the fill tow results in negative 
shear stress between locations F and G, and positive shear between G and H. The maximum 
value of x„ is over twice that of o z in the fill tow. In the warp tows ranges from -0.5 ksi to 
0.5 ksi. In most of the warp tows t„ is very small, which is expected due to the relatively small 
shear modulus (G„) of the warp tows. 

3. 5. 2. 2 R—0.5 Laminate 

Contoui plots of the stress distributions in the undulations of the R=0.5 laminate are 
shown in Figure 56 through Figure 59. The 0 „ stress contour is shown in Figure 56. In con- 
trast to the R =0 laminate (Figure 52), the fill tows in the R-0.5 laminate exhibit both tension 
and compression. The magnitude of o, ranges from a maximum of approximately 95 ksi at 
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Figure 54. Contour plot of in (F/W ? /F) f laminate. R = 0, c 
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Figure 55. Contour plot of in (FfW J /F) | laminate. R = 0, z* = 0.1%. 
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locations A. E, and I. to a minimum of *10 ksi at locations B, D. H. and G. The form of the o, 
distribution in the fill tows is attributed to three factors: first, the axial load of c. - 0.1% con- 
tributes a tensile component of a, to the entire laminate; second, the tendency of the fill tows 
to straighten contributes additional tensile stresses to the fill tows at locations with concave 
curvature (e g. A and I) and contributes compressive components of a. to locations of convex 
curvature (e.g. 0 and G); third, the overall bending of the laminate, as illustrated in 
Figure 50 and Figure 51. contributes additional compressive a, to the fill tows. The net result 
is that the fill tows exhibit areas of tension and compression even though the overall loading 
state is tensile. It is also of interest to note that the maximum value of o, in the R=0.5 lami- 
nate is approximately 95 ksi. which is much larger than the 60 ksi maximum observed in the 
R =0 laminate. Again, this is attributed to the increased bending of the laminate. In the warp 
tows (Figure 56b) o, is less than 1 ksi except in locations A, B. and C. At these locations o„ 
approaches 4 ksi. and is attributed to bending effects. 

The o r stress contour is shown in Figure 57. The maximum value of o r in the curved re- 
gions of the fill tows is approximately 400 nsi. As was the case for the R=Q laminate, local 
maximums and minimums in the o r distributions correspond to their counterpart locations in 
the a, distribution. The maximum value of o r in the R=0.5 laminate (400 psi) is only about 
half that of the R=0 laminate (850 psi). This is attributed to the boundary condition on both 
laminates that F r - 0. In the R=0 laminate the regions of high o„ which generates high o r 
due to Poisson effects, are concentrated a smaller area than in the R=0.5 laminate. 
Therefore, regions of o r are very localized and must be of greater magnitude in order to satisfy 
the boundary condition on F r 

The a z stress contour is shown in Figure 58. The magnitude of o t is approximately 3 ksi 
in the fill tows between locations B and C and between locations J and K. Between locations 
G and H and between I and J, however, a, is near zero. In these regions the upper and lower 
fill tows are bending in the same direction. Thus, the relative z displacement between them 
is small and results in no contribution to the c, component of stress. 
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b) Warp Tows 


Figure 56. Contour plot of O j in (FAW /F) ( laminate. R=0.5, C ] 
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b) Warp Tows 


Figure 57. Contour plot of in (F/W ? /F) ( laminate. R = 0.5, e > =0.1%. 
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The t„ stress contour is shown in Figure 59. It is continuous within the laminate and 
ranges from -15 ksi to about 11 ksi in the fill tows and from -4 ksi to about 2 ksi in the warp 
tows. The distribution in the fill tows is very similar to that of the R =0 laminate. The x„ 
distribution in the warp tows does differ between the R =0 and R=0.5 laminates. In the R=0 
laminate. x B in the warp tows is negligible. In the R=0.5 laminate, however. t„ is nonzero 
due to lack of laminate symmetry. The maximum T n is found in the region between points B 
and C. and results from shear deformation between the fill tows as the laminate deforms. In 
region B-C, t B is approximately -4 ksi. On either side, between A and B and between C and 
0. T n is positive and approaches 2 ksi. 

The stress distributions in the R = 1 and R=4 laminates are illustrated in Figure 107 
through Figure 114. The stress distributions in these laminates are very similar to those 
previously discussed; therefore, they are not discussed in detail. However, the maximum 
stress values found in the undulation are compared for *11 offset ratios in Table 9. As was the 
case for out-of-p ane loading, these results indicate the degree of interaction between the 
undulations. The in-plane normal stress a, appears to be a function of the degree of bending 
of the laminate. The lowest value of a, (69.8 ksi) is found in the case R=0. which exhibits the 
least amount of fill tow bending because the laminate is symmetric about the the x-y plane. 
The maximum value of o, (94.7 ksi) is found in the R=0.5 case, which exhibits the greatest 
amount of fill tow bending of the cases studied. The in-plane normal stress a r and the 
interlaminar normal stress a t are maximum in the R=0 case. The maximum value of o z is 
located in the fill tow at the midplane of the laminate (location G in Figure 54). Location G is 
the point where both of the fill tows in the R=0 laminate will tend towards maximum defor- 
mation in the out-of-plane direction under in-plane loading. Since the fill tow above the 
midplane will tends toward displacement in the +z direction, and the fill tow below fie 
midplane tends toward displacement in the -z direction, a condition of maximum o z develops. 
The distribution of cs y is a result of Poisson effects due to a, and a,. Tne interlaminar shear 
stress x„ does not appear to be a simple function of R. The maximum value (14.3 ksi) occurs 
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A = 0.5 ksi 
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a) Fill Tows 



A = 0 5 Ksi 


b) Warp Tows 


Figure 58. Contour plot of a t in (F/W ? /F) ( laminate. R=0.5, E < = 0.1%. 
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in the R = 1 case and the minimum value (10.2 ksi) occurs in the R=4 case. The magnitudes 
of for the R=0 and 0.5 cases fall between those listed above. 


3.5.3 Interfacial Stress Distributions: (F/W 2 /F), Laminates 


In this section and the following sections on stress distributions in the (FIW^F) w and 
(F/W) z laminates, the stress distributions along the interface of the laminates are discussed 
in detail. As was discussed previously, each ply of the laminate is reinforced in the thru- 
thickness direction due to the woven nature of the fabric; however, no thru-thickness re- 
inforcement exists between adjacent plies. Therefore, the stress state, particularly the 
interlaminar stress state, at the interface betweer adjacent plies is of interest. 

3.5 .3.1 R=0 Laminate 

In-plane Normal Stresses and o y : The interfacial o, distribution in the R=0 laminate is 
shown in Figure 60a. In the far-field region o, is relatively constant at about 0.75 ksi. Small 
increases in a, occur at points B and D, which are within the region of undulation. The dis- 
continuous jump in a, at C can be explained by considering the mismatch in axial modulus. 
E, between the warp and fill tows. Examination of Table 3 shows that £, of the tow used in this 
study is 40.7 MSI and E z is 0.75 MSI. In the far-field region the warp tow is oriented such that 
its axial stiffness corresponds to £* of the tow (0.75 MSI). At point C, however, a„ is evaluated 
in a fill tow, which is oriented such that its axial stiffness corresponds to £, of the tow (40.7 
MSI). Thus, under the in-plane strain loading c, = 0.1% the far-field value of a, in the warp 
tows should be 0.75 ksi and o, at point C should be 40.7 ksi. As indicated previously, the far- 
field value of a , in the warp tows is 0.75 ksi. At point C, however, a, is approximately 38 ksi, 
which is lower than the 40.7 ksi value predicted by CLPT. This difference is attributed to the 
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bending of the fill tows under axial load, which induces compressive or, at C. A contributing 
factor may be the location of the Gauss point at which the stress is evaluated, which is slightly 
offset from C and is at a location where there is slight curvature of the fill tow. As mentioned 
in the discussion of laminate stiffneoC (section 5.3.2), small changes in fiber orientation result 
in significant reductions in axial modulus. The slight increases in o, at points B and D are 
attributed to two factors. First, the reduction in axial stiffness at the undulation due to the 
curvature of the fill tows. At points B and D the contribution of the fill tows to the axial stiffness 
of the laminate is at a minimum due to their orientation with respect to the load axis; therefore, 
the warp tows must absorb a greater share of the load. Second, the deformation of the fill 
tows (Figure 51b) towards the midplane induces compressive a, in the warp tows near the 
undulation. The compressive o t causes a local increase in s x due to Poisson effects and hence 
increases o„. 

The interfacial o y distribution is shown in Figure 60b. As was the case for out-of-plane 
loading, o y stresses are a result of the mismatch in Poisson's ratio between the warp and fill 
tows. In the far-field region o y is evaluated in a We>rp tow and is compressive at approximately 
-100 psi. In the region of undulation o y exhibits a maximum compressive stress of -400 psi at 
points B and D. but becomes tensile at point C where o y is evaluated in a fill tow. The mag- 
nitude of o y at C is approximately 800 psi. The behavior of o y at points B. C, and D is attributed 
to Poisson effects from o , and c x . 

Interlaminar Stresses a ■ and t : The distribution of the interlaminar normal stress o z along 
the interface is illustrated in Figure 61a. The far-field value of o z is zero, which is consistent 
with classical laminated plate theory (CLPT). At points B and D, o z is compressive (approxi- 
mately -2.5 ksi). At point C, <j z is tensile (approximately 5 ksi). The o z distribution at the un- 
dulation is explained by considering the deformation of the fill tow (Figure 51). As discussed 
previously, the tendency of the fill tow is to straighten under axial load. This induces tensile 
stresses at C and compressive stresses at points B and D.. 
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The interlaminar shear stress x„ (Figure 61b) is essentially zero throughout, which is 
consistent with CLPT. The exception occurs on either side of C. where x„ exhibits values of 
0.3 and -0.3 ksi. This is a point of stress concentration and is a result of the dealized geom- 
etry of the model. The oscillation of about x„ = 0 in the region of undulation is attributed 
to numerical instabilities associated with the finite element method. These instabilities are 
often associated with regions of stress singularities, where linear elasticity predicts infinite 
stresses. 

3 . 5 . 3.2 R=0. 5 Laminate 

In-plane Normal Stresses o % and o y : The distribution of the in-plane normal stresses o x and 
o y along the interface for the R=0.5 laminate are illustrated in Figure 62. The far-field value 
of a, is approximately 0.75 ksi, which is consistent with CLPT. The distribution of g„ in the 
region of undulation is relatively complex. At points C and E, g, is evaluated in fill tows and 
is approximately 3.5 ksi. The increase in o, at points B and F is attributed to the deformation 
of the fill tows (Figure 51). As was the case for the R = 0 laminate, the tendency of the fill tows 
to straighten under axial load results in compressive g z at locations B and F. This induces 
Poisson expansion of the warp tows in the x direction and thus increases o„ at these points. 
The distribution of o x at location D is attributed to bending of the laminate. This laminate is 
unsymmetric about the x-y plane, with maximum asymmetry existing at location D. Therefore, 
bending of the laminate is greatest at D and a, due to bending is maximum. 

The Gy distribution in the R=0.5 laminate is similar to that of the R=0 laminate. The 
far-field region is m compression and is approximately -125 psi. At locations B and F, o y is 
approximately -300 psi. As was the case for the R = 0 laminate, this behavior is due to defor- 
mation of the fill tows, which increases o y at B and F At locations C and E, G y is evaluated in 
fill tows and is approximately 300 psi. 
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Figure 60. Interfacial ct > and O y In (F/W 2 /F) ( laminate. R = 0, = 0.1%. 
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Interlaminar Stresses o and x : The distributions of the interlaminar stresses a. and t„ are 
illustrated in Figure 63. In the region of undulation a z is tensile at locations C and E and 
compressive at points B, D, and F. As was the case for the R=0 lammate, this behavior is 
attributed to the tendency of the fill tows to straighten during deformation, which induces 
tensile stresses at the peaks of the undulation (C and E) and compressive stresses away from 
the peaks (B, 0, and F). 

The interlaminar shear stress is nonzero throughout the undulation region. Unlike the 
R=0 laminate, the R=0.5 laminate is unsymmetric about the x-y plane. Therefore, bending 
of th*' laminate occurs under axial strain and shear stresses result. The maximum is -4 
ksi and occurs at location D, which is midway between the undulations. Point 0 is the location 
of maximum asymmetry; therefore, it is the location of maximum shear stress. Positive values 
of (= 2 ksi) occur midway between points B and C, and between points E and F.. This 
behavior is also due to bending and shear deformation of the laminate. 

3. 5. 3. 3 R~1 and 4 Laminates 

The interfacial stress distributions for the R — 1 and R~4 laminates are presented in Ap- 
pendix A (Figure 125 through Figure 128). The form of the stress distributions in these lami- 
nates are similar to those of the R = 0 and R = 0.5 laminates. Similarities and differences 
between the laminates are summarized in section 3.7. 

3.5.4 Interfacial Stress Distributions: (W/W,/F) w Laminates 


The interfacial stress distributions for the (F/Wj/F),,, for R = 0 and 0.5 laminates, under 
in-plane loading are shown in Figure 64 through Figure 66. As previously mentioned in the 
discussion on out-of-plane loading, these laminates differ from the (F/I/Vy^, laminates in that 
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the warp tows in the (FfWJF)., laminates exhibit curvature and bound the midplane, whereas 
in the (FIWJF), laminates the fill tows exhibit curvature and bound the midplane (Figure 18). 


3 5.4.1 R~0 Laminate 

In-plane Normal Stresses 0 ( and o*: The a, and <J r stress distributions in the R =0 laminate 
are shown in Figure 64. In the far-fielii region <j, is evaluated in a warp tow, which exhibits 
an axial modulus £, of 40.7 MSI. The magnitude of o, is approximately 39 ksi. which is less 
than the 40.7 ksi predicted for an in-plane strain c, - 0.1% using CLPT. This difference is due 
to the presence of the undulation, which is a region of low axial stiffness in the laminate. 
Therefore, there is nonuniform z, distribution along the laminate, with the region of undulation 
under an axial strain c, > 0.1% and the far-field region under a strain z, < 0.1% . In the re- 
gion of undulation, a, is evaluated in a fill tow, which has a low axial modulus (0.75 MSI). The 
magnitude of <7, in the undulation is approximately 0.75 ksi. 

The a, component of stress is tensile and approximately 100 psi in the far-field region. 
At location C, o r is compressive and approximately -600 psi. At locations 8 and D. o r is ap- 
proximately 700 psi. The o, at B and D is attributed to large o r -induced Poisson extension 
of the warp tows in the y direction. 

Interlaminar Stresses o and t : The distribution of the interlaminar normal stress cr, is shewn 
in Figure 65a. At point C. n z is -4 ksi. The laminate is in compression at C because the 
tendency of the warp tows to straighten under load causes the peaks of the undulation to 
move towards the midplane, putting the central warp region under compression. At locations 
B and D the tendency is for the warp tows to separate, resulting in tensile o r . 

The interlaminar shear stress t„ {Figure 65b)is zero in far-field region, but is nonzero in 
the region of undulation. This laminate is symmetric about the x-y plane and thus should ex- 
hibit zero i„ everywhere along the midplane. As mentioned in the discussion on the 
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[F/Wj/F)„ R=0 laminate, this discrepancy is attributed to the presence of stress singularities 
(points B and 0 in the ( F/WJF ).. R =0 laminate) and the numerical instabilities associated with 
them. 


3JS.4JI R=0.5 Laminate 

In-plane Normal Stresses o f and o y : The interfadal o, distribution in the R=0.5 laminate is 
shown in Figure 66a. In the region of undulation the a, distribution is complex, and is best 
understood by considering the deformation of the laminate (Figure 67). Figure 67 indicates 
that the laminate exhibits severe bending, which results in a radically different a, above and 
below the interface. Below the interface, o, is high at B and above the interface o, is high at 
F (Figure 66a). The low value of a, at location O is due to the change in material properties. 
At point D. o, is evaluated in a fill tow. which has a low axial modulus. 

The interfacial o r distribution (Figure 66b) follows that of the o, and <j 4 distributions due 
to Poisson effects. As was the case for a, the c r distribution is different above and below the 
interface. In the far-field region o, is tensile and approximately 100 psi. At locations C and 
E, which are above and below the interface, respectively, c r is compressive and is approxi- 
mately -250 psi. At locations B and F. u r is tensile and approximately 350 psi. 

Interlaminar Stresses and t u : The interlaminar normal stress o z is shown in Figure 68a. 
As in the laminates discussed previously, the a, distribution can be related to the tendency for 
the warp tows to straighten under axial strain. In this laminate the deformation of the warp 
tows results in compressive o z at location D and tensile o t at locations B. C. E, and F. The 
maximum tensile o, is approximately 1.8 ksi. at points B and F. The maximum compressive 
o, is approximately -1.2 ksi. at point D.. Continuity of a, does not appear to be satisfied at 
locations B, C. E. and F. Locations B and F are singular points, which are not accurately 
modeled by the finite element method. Locations C and E are regions of high stress gradients. 
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Figure 65. Interfacial O* and In (F/W 2 /F)^ laminate. R=0, z* = 0.1%. 
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which require a denser mesh than that used in this study in order to be accurately predicted. 
The interlaminar shear stress z a (Figure 68b) is induced by the bending of the laminate under 
axial load, and ranges from +4 ksi at location D to -5.5 ksi at locations D and F. The x„ dis- 
tribution indicates that continuity is satisfied at the interface. 

3. 5. 4. 3 R—1 and 4 Laminates 

The interfacial stress distributions for the R = 1 and 4 laminates exhibit many of the same 
characteristics those already discussed, and the physical arguments used to explain their 
behavior are similar. As such, further discussion is left to a later section where the similari- 
ties and differences between the laminates are summarized. The distributions, however, are 
presented in Appendix A. 


3.5.5 Interfacial Stress Distributions; (i-7W) 3 Laminates 


As may be recalled from the discussion of the out-of-plane loading results, the stacking 
sequence of the (F/W) t laminate differs from that of the (F/V^/F) laminate in that the ( F/W) z 
laminate is not symmetric in the far-field region. This results in the midplane being bounded 
by f jws of opposite orientation. The tow above the midplane is a warp and the tow below the 
midplane is a fill. 

3.5.5.1 R—0 Laminate 

In-plane Normal Stresses and a : The interfacial o, distribution in the R=0 laminate is 
shown in Figure 69a. The far-field value of a, in the warp tow (above the interface) is ap- 
proximately 0.75 ksi, and agrees with CLPT. In the fill tow (below the interface) a, is approxi- 
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Figure 67. Deformation of (F/Wj/F)^, R = 0.5 laminate, t^-0.1% 
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mately 34 ksi, which is far less than the 40.7 ksi predicted by CLPT. As was the case for the 
(F/Wj/F)* laminate, this result is due to the presence of the undulation, which is a region of low 
axial stiffness and hence causes a nonuniform e„ distribution to exist in the laminate. Conse- 
quently. e, < 0.1% and o„ in the fill tow is reduced. The effect is negligible in the warp tow 
due to its low axial modulus. The o, distribution in the region of undulation is more complex. 
At point C, e, is compressive at about -4 ksi. This is due to bending of the laminate, which 
leaves the fill tow in compression at C. The slight increase in o, at locations B and D is also 
attributed to bending of the laminate, which increases o, in the warp tows near the undulation. 

The o y distribution {Figure 69b) indicates that in the far-field regions the warp tow (above 
the interface) is in tension and the fill tow in compression. The magnitudes of a y in the far-field 
regions of the warp and fill tows are 100 psi and -125 psi, respectively. This indicates an ap- 
parent lack of force equilibrium in the far-fteld region. However, in an unsymmetric laminate 
neither stresses nor strains are constant within each tow; therefore, midplane stresses are 
not an accurate indicator of force equilibrium in the laminate. The a y distribution in the region 
of undulation can be related by Poisson effects to the a* and o z distributions. 

Interlaminar Stresses ct j and The interlaminar normal stress o z (Figure 70a) is maximum 
and tensile ( = 2 ksi) at location C. At locations B and D, ct z is compressive at approximately 
-1.2 ksi. As was the case for the (F/Wj/F)„ and ( F1WJF ), laminates, the a z distribution is due 
to the tendency of the fill tows to straighten under the applied strain. This induces tension at 
the peak of the undulation (C) and compression c i either side (B and D). The far-field value 
of t„ is zero, which is consistent with CLPT. 

The interlaminar shear stress (Figure 70b) is 4 ksi and -4 ksi at locations B and D, 
respectively The shear stress is a r»iult of shear deformation which in turn results from the 
asymmetry of the laminate in the region of undulation. The far-field value of t„ is zero, which 
is consistent with CLPT. 
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3.5 .5.2 R—0.5, 1, and 4 


The stress distributions in the R-0.5. 1, and 4 laminates can be described using the 
same physical arguments used to describe the laminates discussed previously; therefore, the 
specifics of the distribution of each component of stress for each of the remaining offset ratios 
will not be presented. However, the interfacial stress distributions for these laminates are 
presented in Appendix A. A summary of the results for in-plane loading is presented in Sec- 
tion 3.7. 


3.5.6 Influence of Undulation Aspect Ratio n 


The influence of the undulation aspect ratio £2 on the maximum interfacial stress was 
evaluated for a (F/W/F), laminate with a.' 1 offset ratio R = 0.5. The maximum values of o , , o y , 
o z , and t c along the interface, as a function of £2, are presented in Figure through 
F.^.re 74. The variation of a, (Figure 71) can be divided into two distinct re b ,«.,,6. For undu- 
lation aspect ratios £2^ 12. the maximum o, occurs at location C and increases approximately 
linearly with increasing £1 This can be explained by recalling from the discussion of the out- 
of-plane loading results that as £2 increases the geometry of the laminate approaches that of 
a symmetric cross-ply laminate, which does not bend under axial load. Therefore, a, will in- 
crease since no compressive a, due to bending exists at C. At Si values less than 12, the 
curvature of the fill tow is large enough that significant bending occurs under applied axial 
strain and location C is in compression; the location of the maximum tensile a, changes from 
location C to location D. The o y distribution (Figure 72) shows that o y decreases with in- 
creasing Si. At £2 = 6, o y is approximately 0.46 ksi and decreases to 0.22 ksi at £2 = 24. The 
form of the o y distribution follows that of the o, distribution, which is shown in Figure 73. 

The interlaminar normal stress a, (Figure 73) is maximum ( = 3.1 ksi) at £2 = 6 and de- 
creases to ipproximately 1 ksi at £2 =24. This is also due to the reduction in curvature of the 
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Figure 69. Intcrfacial o > and in (F/W) 2 laminate. R = 0, e =0.1%. 
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fill tows with increasing £2. As £2 increases the response of the hminate approaches CLPT, 
for which o x is zero. The interlaminar shear stress t a (Figure 74) does not vary significantly 
with £2. The magnitudes of x a ranges from 3.9 ksi at ft =9 to 3.4 ksi at £2=24. which is a dif- 
ference of approximately 15 percent. The form of the distribution is attributed to the ori- 
entation of the curved portion of the fill tows with respect to the load axis. At £2 = 9 the fill tows 
at the undulation are oriented approximately 45° with respect to the load axis, which is the 
condition for maximum shear. Thus, at £2^ 9 the magnitude of x a is expected to decrease. 


3.6 Results: Thermal Loading 


3 R.1 Stress Distributions In Region Of Undulation 


The thermal stress distributions (per °F temperature change) in the undulations of a 
(F/Wj/F),, R =0.5 laminate are presented as contour plots in Figure 75 through Figure 79. The 
in-plane normal stress o„ is tensile and approximately 10 psi in the far-field regions (D and 
E) of the fill tow (Figure 75a). The far-fieid a, in the warp region (Figure 75b) is approximately 
-10 psi, indicating satisfaction of force equilibrium in the x direction. These stresses are due 
to the mismatch in the coefficients of thermal expansion between the .'ill tows and the warp 
tows. The CTE of the fill tows is much less thar that of the warp tows (Table 3); therefore, 
under a positive AT the fi'l tows are in tension and the warp tows in compression. The o„ 
stress contour in the region of undubtion is more comple/ and is attributed to the deforma- 
tion of the Til 1 lows ur.der thermal load, which is analogous to the straightening of a curved 
bear. (Figure 76) Thu response induces tension in the fill tows at locations A and H (25 psi) 
and axial compression in the regions between B and C and between F and G, where o, ranges 
horn -5 to -15 psi The warp region is in axial compression throughout. 
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The o y contour (Figure 77) shows that o y varies little within the undulation. In the fill tows 
o y is compressive and ranges from -9 to -10 psi. In the warp region a y is tensile and ranges 
from 9 to 10 psi. satisfying equilibrium. This result is also due to the mismatch in CTE's be- 
tween the warp and fill tows. In the y direction the CTE of the fill tows is greater than that of 
the warp tows. Therefore, under uniform positive AT the fill tows are under compression in 
the y direction and the warp tows are in tension. 

The contour of the interlaminai normal stress o z is shown in Figure 78. In the far-field 
regions a z is zero for both the warp and fill regions, which is consistent with CLPT. A large 
gradie.it in o z exists in the regions B-D and E-G. Maximum tensile o z is found in the regions 
between C and D and between E and F, where o t is approximately 1.5 psi. These regions 
correspond to the regions in the deformed laminate where the fill tows are most oriented to- 
wards the z direction (C-D and E-F in Figure 76). Therefore, the stiffness in the z direction is 
greatest and hence o z is maximum. At locations B and G, a z is compressive (approximately 
-3 psi). This result is attributed to bending of the fill tows under thermal load, which forces the 
fill to towards the midplane of the laminate at these points, and results in compressive o z . 
Inspection of the o z distribution indicates that continuity is satisfied. 

The contour of the interlaminar shear stress x^ is shown in Figure 79. In the far-field 
regions (A and F) x„ is zero as expected. Within the undulation, however, r n is nonzero due 
to the asymmetry of the laminate. As was the case for in-plane loading, lack of symmetry in 
the laminate induces shear deformation in the undulation and hence shear stresses. In the 
fill tows (Figure 79a) x a is negative (-3 to -P psi) in the regions B-C and H-l and is positive (1 
to 3 psi) in the regions C-E and G-H. The shear stress in the warp tows (Figure 79b) is positive 
in regions B-C and D-E at about 1 psi. Between the undulations (region C-D) x„ is negative 
and ranges from about -2 psi to -3 psi. 
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3.6.2 Interfacial Stress Distributions 


3.6.2.1 In-Plane Normal Stresses a and a 

* » 

The interfacial thermal stress distributions are presented in Figure 80 and Figure 81. 
Again, these stresses result from the mismatch in thermal expansion coefficients between the 
warp and fill tows. The in-plane normal stress a, is compressive ( = —9.5 psi) in the far-field 
region. At locations C and E, a, is evaluated in fill tows and is tensile at approximately 12 psi. 
The o y distribution (Figure 80b) shows that the warp tows (far-field region) are in tension (9.5 
psi) and the fill tows arc in compression ( -9.5 psi at C and E). This indicates satisfaction of 
force equilibrium in the y direction. 

3. 6.2.2 Interlaminar Stresses o and x 

* » 

The interlaminar normal stress cs t distribution (Figure 81a) indicates that o z is tensile at 
the undulation peaks (C and E) and compressive on either side of the undulations (B. D.and 
F). This distribution is attributed to the deformation of the laminate under thermal load. The 
deformation of the fill tows is such that they tend to separate from the midplane at locations 
C and E compress the midplane at locations B, D and F. The interlaminar shear stress, x a , 
distribution is presented in Figure 81b. Positive shear stresses of approximately 1.6 psi exist 
at locations B and F. At location i), which is the location of maximum asymmetry, is ap- 
proximately -3.1 psi. In the far-field region x„ is zero, which is consistent with CLPT. Both the 
a, and x„ distributions indicate that stress continuity is satisfied at the interface. 
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a) Fill T ows 



A-0.5p3i/°F 

b) Warp Tows 


Figure 75. Contour plot of 0^ in (F/W } /F) f laminate under thermal load. R = 0.5, AT = + 1°F. 
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Undeformed 

Deformed 


Figure 76. Deformation of fill tows in (F/W ? /F) f laminate under thermal load. R=0.5, AT = +1°F. 
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3.7 Summary and Discussion 


3.7.1 Stiffness 


Woven fabric composites are unique in that laminate stiffness varies along the axial co- 
ordinate. The far-field region, which has the geometry of a cross-ply lamir.ate. exhibits elastic 
moduli which can *■ a predicted from classical laminated plate theory (CLPT). In the regions 
of undulation the moduli are a function ' ' .tacking sequence and offset ratio R (Table 5). The 
in-plane modulus E. is most influenced by the above geometry, parameters, with reductions 
of greater than 5C% of the far-field modulus predicted for the cases considered. The moduli 
£, and E z are relatively independent of stacking sequence and offset ratio. 

The effect of the undulation on £„ of the unit cell is presented >n Table 6. The results ln- 
oicate that undulations reduce the axial modulus. £„ up to 18% from that of the far-field value. 
The average reduction in £, is about 12%. and is relatively independent of R. 


3 7.2 Out-of-Plane Loading 

3.7.2.1 Region of Undulation 

The results of the out-of-plane loading analysis show that the stress distribution in the 
laminates s relatively complex, and can be related to the two distinct regions of interest that 
exist in each stacking sequence ((F/W r 'F) t ( F/W r 'F) w . and (F/W),) for all offset ratios R In the 
far-fielcf region the geometry of the laminate is analogous to a cross-ply laminate. As such, 
the normal components of stress o„ o,. and o t are constant throughout each tov: and the 
interlaminar shear stress r, z is zero. In the region of undulation, all components of stress are 
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nonzero and are a function of stacking sequence, offset ratio R. and undulation aspect ratio 
£1. The stress distributions in the region of undulation can be explained in terms of local de- 
formation of the tows, which is analogous to bending and shear of a beam (rig 24). Under 
out-of-plane loading the bending of the tows generates large a, as well as increases the 
magnitudes of a r o r and x n c *er the far-field values. 

2.722 Interfacial Stresses 

The stress distribution along the interface between the two woven fabric lamina is of 
special interest because interlaminar failure was observed experimentally (Chapter 2). The 
maximum interfac'al stresses, as a function of stacking sequence and offset ratio R. are 
summarized in Figure 82 through Figure 85. 

In-plane Normal Stres. es o* and o^: The in-plane normal stress o, (Figure 82) is largest in the 
{F/Wf/F), laminate and smallest in the (F/Wj/F)_ laminate. This can be explained by consid- 
ering the geometry of the laminates as illustrated in Figure 14. In the {F/WJF), laminate the 
convex side of the fill tow. which exhibits maximum o,. lies along the interface In the 
{F/Wj/Fj m laminate, however, the convex side of the tow lie along the top and bottom surfaces 
of the laminate. Thus, in the {F/W t 'F) w laminate the region of maximum a, is not part of the 
mtertacial stress distribution. The o, in the (F/VV)j laminate differs from the other laminates 
in that the maximum value of o, is found in the R =0 5 case This is attributed to the geometry 
of that particular laminate Apparently, the orientation of the fill tows is such that naximum 
interaction, in the form of bending between the two undulations, occurs, and results in high 
o.. It should also be noted that o, varies little with offset ratio for R i 1 

The maximum a, (Figure 83) is found in the [F/W l /F) m laminate, and is approximately 10% 
higher than the largest o r predicted for both the (F1WJF), and [F/W) z laminates Additionally. 
o y is essentially independent of R. 
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Interlaminar Stresses o and r : The variation of '-e interlaminar normal stress o. 

I at * 

(Figure 84) exhibits the same characteristics as the o r variation. The maximum o, is found in 
the (F/MT,/F). laminate and is nearly independent of R. It should also be note ' that o, is rela- 
tively independent of axial coordinate. The maximum o t predicted for the region of undulation 
was less than 15% greater than that predicted for the far-held region. 

The variation in the interlaminar shear stress (Figure 85) differs from the other com- 
ponents of stress in that t„ is strongly fien»n riant on laminate stacking sequence. The largest 
t,. is found in the (F/My'F), laminate and the smallest is found in the (F/W/F), laminate. 
For ail three cases studied is essentially independent of R for R ^ 1. The variation in t„ 
for R < 1 is attributed to lelative changes in shear deformation of the laminate associated 
with the proximity of the undufadons to each other. 

Recognizing that the (F/W/F), and (F/wyF)_ stacking sequences represent the two 
-ross-sections of the (FW^F) laminate (Figure 14). it appears that the (F/W), laminate would 
be least susceptible to interlaminar failure due to the low predicted. However, short-beam 
shear and transverse tensile strength tests of several unidirectional carbon-carbon materials 
conducted by Perry and Adams 4 indicate shear strengths approximately five times greater 
than the transverse tensile strength Consequently, it is predicted that a. . which is not a 
strong function of laminate geometry, dominates failure, hence, all the stacking sequences 
evaluated are equally susceptible to interlaminar failure. 


3.7.3 In-Plane Loading 


3.7. 3.1 Region of Undulation 


Similar to the out-of-nlane loading case, the stresses due to in-plane loading can be de- 
scribed in terms of the geometrv of the laminate, namely the region of undulation and the 
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Figure 84. Maximum interfaciai O j as a function of offset ratio and stacking sequence. r, =0.1%. 
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far-field region. In the far-field region the interlaminar stresses o z and x„ are zero and the 
in-plane normal stresses o, and a r are essentially those predicted from CLPT. In the region 
of undulation, however, all components of stress are nonzero and. like the out-of-plane loading 
case, can be related to the bending of the tows. Under in-plane loading the tendency is for 
the tews to straighten (Fig. 68). Thus, on the convex side of the tows o, < 0 and on the con- 
cave side a, > 0. Additionally, significant o, and are generated due to the orientation of 
the curved tows in the undulation, which is analogous to an off-axis unidirectional laminate 
under tensile load. 

3.7.3.2 Interface I Stresses 

The maximum interfacial stresses as functions of laminate and offset ratio R are sum- 
marized i:i Figure 86 through Figure 89. 

In-plane Normal Stresses ind cy The in-picne normal stress a, (Figure 86) is largest in the 
(F/W*/F) w laminate and smallest in the (F/WJF), laminate. Also, a, is relatively independent 
of R for R i 0.5 for ail three laminates. The radical difference in a, for R=0 is attributed to 
the degree of symmetry of the laminates. The {F/WJF),, and (F/W^F), laminates are symmetric 
about the x y plane Thus, no bending of the laminate occurs during in-plane loading For the 
(F'M/j/F)„ laminate the bending contribution to a, . which is tensile, is eliminated. In the 
(F/Wj/F), laminate the curved tows are oriented such that under in-plane load compressive 
o, is exhibited at the interface between the plies. Therefore, for the case R = 0, compressive 
o, is eliminated due to symmetry. In the (F/W) 2 laminate the tows nest such that the laminate 
exhibits the lowest axial modulus E, of those analyzed; consequently, for an applied axial 
strain c„ a, is significantly reduced. 


Finite Element Stress Analysis 


155 



The a r distribution (Figure 87) shows that for all three laminates, o y is relatively inde- 
pendent of R i 1. For R < 1, a y follows the trend of o, for the (FIWJF), and (F/W) 2 laminates. 
For the (F/wyF) w laminate, follows the trend of the o* distribution. 

Interlaminar Stresses a and x : The variation of the interlaminar stress a., as a function of 

I xz 1 

stacking sequence and offset ratio, is presented in Figure 88. In the {F/WJF), and (F/Wj/F),, 
laminates the maximum o z is found at R=0. In the (F/W) 2 laminate a, is minimum at R=0. 
For R s 1, o, is relatively independent of R. It should also be’ noted that the maximum 
interfacial o z generated by the applied in-plane strain z, = 0.1% is much greater than o t 
generated by the applied out-of-piane strain t, = 0.1% . Under in-plane loading the maximum 
c z is approximately 5.1 ksi; under out-of-plane loading the maximum o z is approximately 0.9 
ksi (Figure 84). 

The maximum interlaminar shear stress x„ (Figure 89) is found in the (F/VVj/F), laminate. 
It ranges from approximately 0.4 ksi at R = 0 to about 6.8 ksi at R = 1. In the (F/Wj/F),, laminate 
T„ ranges from aboui 0.2 ksi at R = 0 to approximately 4 ksi at R = 0.5. In the ( F!W) Z laminate, 
x,j ranges from 4 to approximately 5 ksi. As was the case for the other components of stress. 
x„ is essentially independent of R for R £ 1. 

The magnitudes of the interlaminar shear stresses indicate that the regions of undu- 
lations are possible sites for delamination. The interlaminar normal stress, o r , which ranges 
from 2 ksi to 5 ksi for the cases studied, is much larger than the 1 ksi out-of-plane tensile 
strength reported by Maahs **. Additionally, the predicted x„ is large, ranging from 2 ksi to 
almost 7 ksi. This is consistent with the results ofWagnecz, who observed delamination at the 
undulation in an incrementally loaded Celion 3000/PMR-15 graphite/polyimide woven fabric 
composite 3 ' 
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3.7.4 Thermal Loading 


The analysis of the (F/WJF),, R =0.5 laminate under thermal load conditions indicates that 
significant thermal stresses can be generated during processing and application. The mag- 
nitude of the predicted a, is on the order of 25 psi/'T. fhe magnitude of the predicted o y is 
on the order of 10 psi/°F. The interlaminar thermal stresses a z and x„ generated are on the 
order of 3 psi/°F and 6 psi/°F, respectively. 

The possible failure modes associated with thermal loading may differ depending on 
whether the laminate is undergoing positive AT (e.g. an increase from room temperature to 
application temperature) or negative AT (e.g. a decrease from fabrication temperature to room 
temperature). A positive AT will result in the stress distributions exhibited in Figure 80 and 
Figure 81 such that both the warp and fill tows will be in axial tension (o, > 0) and transverse 
compression (o 2 < 0). Consequently, the direction of maximum tensile stress will correspond 
to the direction of maximum strength and the direction of minimum strength will be in com- 
pression. Failure, should it occur, would be expected to be due to axial breakage of the tows. 
The distribution of the interlaminar normal stress a z (Figure 81) indicates that under positive 
AT delamination may occur at peaks of the undulations (locations C and E). The nonzero 
interlaminar shear stress x« can also be expected to contribute to delamination. If the lami- 
nate is subjected to negative AT the interfacial stress distributions will be of opposite sign to 
those shown in Figure 80 and Figure 81. Thus, the sign of the in-plane normal stresses will 
be such that the tows will be in axial compression and transverse tension. This indicates that 
transverse cracking of the tows can be expected, which is consistant with the observed 
microstructure of the material used in the out-of-plane tensile tests (Figure 5). In those 
micrographs the crack density was observed to be approximately 100 cracks/inch. The dis- 
tribution of the interlaminar normal stress a z indicates that o x will be tensile at locations B, 
D, and F. Thus, under negative AT delamination could be expected to initiate at those sites. 
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Figure 87. Maximum interfacial as a function of offset ratio and stacking sequence, r. 


= 0.1%. 


59 









4.0 Effect of Fiber Anisotropy on Thermal Stresses 
in Fibrous Composites 


4.1 Introduction 


The extreme thermal environment seen by carbon-carbon composites during fabrication 
and am .cat'O-n requires more careful consideration of the choice of fiber than would be re- 
qu red for apmtes.-poxy cc.npo'- tes used in room-temperature applications. In most 
graphife-epox> composite applications T he choice of fiber is based on fiber strength and 
tfnes: in the ax. al direction h applications where t..ermal loads exist the a v ial coetticient 
cf the.c.ial expansion also becomes important In carbon-carbon composites, however, nigh 
thermal loads and t.ie presence of a brittle matrix suggests th«.t. in add.tion to axial properties, 
the transverse fiber properties can significant , influence the stress distribution in the com- 
j •• .r and hence greatly mtluence strength This chapter addresses the issue of tne micro- 
structure of th-> f 1 o ' .Is influence on transverse elastic constants, and the lesulting variations 
in tnermal stress distributions 
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High axial stiffness in graphite Fibers is obtained by processing the Fiber precursor such 
that the stiff basal planes of the graphite crystals Figure 90 are evented nearly parallel to ihe 
longitudinal axis of the fiber. “ * In the transverse direction, however, the orientation 

of the basal planes can result in many different microstructures. Examples of four types of 
microstructures observed in graphite Fibers are shown in Figure 91. In Figure 91a the basal 
planes a<e arranged circumferentially around the Fiber. This structure is commonly called an 
'onionskin" structure, and it would be expected that E, > £ r In Figure 91b the basal planes 
are arranged radially, for which E, > £,. Figures 91c and d show combinations of radial and 
circumferential microstructures. Figure 91c shows a radially oriented core with an onionskin 
st.aath. Figure 91d shows a random core with a radially oriented sheath. The structures 
shown in Figure 91a and Figure 91c are normally associated with polyacronitrile-based (PAN) 
fibers and the structures shown in Figure 91b and Figure 91d are commonly associated with 
pitch-base fibers. 37 :f “ * " 

The elasticity (emulation vsed in this paper follows the previous works of St. Venant 
Voigt Mitinskii *\ Lekhmtskn 45 **, and Cohen ei al. 47 ** 43 Lekhmlskii provides the form of the 
solution for a variety of loading conditions on solid and hollow ivluy^rs possessing 'cylin- 
drical anisotropy". 45 Cohen and co-workers provided explicit forms cf the equations for tsm:- 
nated composite tubes subjected to thermal loading 

Ip this chapter, the anisotropic elasticity solution is used to explore the interrelationship 
between fiber orthotropy. as exhibited by the transverse micr^structures. and the stress dis- 
tributions in a fiber-matrix composite under uniform thermal lead As will be shown, the type 
of orthotropy radically affects the thermal stress distribution in the fiber In particular, radial 
orthotropy (C„ > C*J in the center of the fiber (Figure S 1 b-c) results in singular stresses at the 
center of the f.oer for all three rormrl components of stress This has obvious negative con- 
sequences for th development of damage m the form of fiber splitting and fibc- - breakage 

’’he results of this analysis provide helpful insight into the structural integrity of the fiber 
a, a function of microstructure. Such insight may pro/e helpful in choosing a fiber for ,i par- 
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Figure 90. The Graphite Crystal Structure 
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Figure 91. Transverse Microstructures o( Graphite Fibers 
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ticular application, such as in carbon-carbon composites, which are subjected to large the imal 
loads. In addition, the results may aid in guiding fiber development for improved properties. 


4.2 Mathematical Formulation 


Consider a long fiber of radius a in an isotropic matrix of thickness b - a under uniform 
thermal load (Figure 92). 

Oue to 3xial symmetry the hoop displacements are zero, stresses and strains are inde- 
pendent of 0. and there s no shear-extension coupling. 

Therefore, the thermoeiastic stress-strain relations are 




*-xjc < -'x9 C xr 


f -x ~ a iA T ) 

n 0 


Cix C 03 Cf) r 


H ~ OqAT > 



C rx Crr 

1 

r, r - a# A T ) 


z r0 T xr ®xrYjrr* *Qx ^OxfOx 


[4.1] 


where C, t are stiffness coefficients and [C] represents the appropriate matrix for the three 
normal components of -tress, a are coefficients of thermal expansion, and AT is the uniform 
temperature change. The nonvamshmg equilibrium equations for this axisymmetric problem 
are 


i3< J r 
tr 


+ 


T(°r ~ ° 0 ) 


0 


[4.2a] 


Tr 



The strain-displacement relations can be written 


[4 2b] 
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[43a] 





r -» = f 

[435] 

C = 

[4.3c] 

cx 

Trfl = 0 

[4.3cf] 

^ 

[4.3e] 

cr 

f8x = 0 

[4.3f] 


where u. v. and w are axial, hoop, and radial displacements, respectively. 

Substituting Eqns. 4.3a-f into Eqn. 4.1 and substituting the resulting equations into the 
equilibrium equation 4.2a yields the governing differential equation 


( w 


Cr 2 


+ ± 


cw 

cr 


- Coo-y = j-( c ox ~C rx K + 7“( c r ; - C 0J )aiAT 


[4.4] 


for generalized plane strain with uniform axial strain c,. Here and throughout the paper re- 
peated sub cnpts ij are summed over x. r, and 0. Equation 4 may be solved for the case of 
apolied thermal load or the case of applied axiai strain. 

The general solutions io Eqn. 4 fallowing Cohen and Hyer's 4 * treatment of an orthotropic 
tube under umlorm thermal load are: 


a. transversely orthotropic fiber (C.„ * C,.), 

w(r) = A^' + A 2 r ^ + H^i: x r -t H 2 ATr [4 5a] 


where 
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Hi = 


Q)X CfX 
C rr ~ Coo 


[456] 


_ (C f , Ca,)a ( 

W 2 


(C rr c 00 ) 


[4.5c] 


b. transversely isotropic fiber (C„ = C„), 


w(r) = 4,^ + + G 1 e x t Inr + G 2 A7r In r 


[4.6a] 


where 


G, = 


Q)x ^rx 


2 C, 


00 


[4.66] 


q _ (Cri C 0l )a, 

2 2C eo 


[4.6c] 


In both of the above displacement fields (4.5a and 4.6a), X 1? are defined 


- 1,2 


= ± 


-On 


10.5 


[47] 


The solution of the second equilibrium equation (4.2b) is 




[48] 


where K is a constant of integration determined from the boundary condition on t„ at r = 6 
For a composite with an orthotropic fiber and an isotropic matrix the equations for the 
radial displacements are 


w f (r) = A{/' f A 2 r X * + h[t./ + H 2 ATr 


[4 9a] 
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w m (r) = A™r + A 2 ~ 


[4.96] 


where the superscripts fandm refer to the fiber and matrix, respectively. = ± 1. and 
Gr = Gf = 0 for a material which is transversely isotropic in both elastic and thermal con- 
stants. 

The equations for the normal components of stress, for the fiber and matrix, can be ob- 
tained by substituting Eqns. 4.5-4 6 into Eqns. 4.3a-b and then substituting the results into the 
constitutive relations (Eqn. 4.1). The resulting equations are: 


a. transversely isotropic fiber ( C' rr = C' w , A u = ± 1), 

o' = A f ,(c!v + C! r ) + A^ - Cl) ±. + C^ x - ctfAT [4. 10a] 

r 


b. transversely orthotropic fiber ( C'„ * CJ,). 


o' = A' (Cl + C^.'jr 1 ' - 1 + A[(cl + Cy 2 )/-* + l\z x + A/,'A7- [4.10b] 


where 


L' = Cl + H\(CI + Cl) 

n!=H^(cI + cl) - c'af 


c. isotropic matrix, 


m _ .m, m . A m, r m -rru 1 

O, ~ (C,o + C, r ) + A 2 (C l0 ~C ir ) — 


[4.10c] 


+* C ix r. x C t j (if AT 


The five constants A [, Al. A?, A?, andi;„ are determined from the following five conditions. 
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1. The radial displacement iv must be zero at r = 0. This condition, plus the fact that 
XJ < 0 (Eqn. 4.7), requires A£ be zero to avoid a singularity in w at r = 0 for both 
transversely orthotropic and transversely isotropic fibers (Eqn. 5.5 and 5.6). 

2. Continuity of w at the fiber-matrix interface w'(a) = w m (a) requires that: 

a. for a transversely orthotropic fiber ( C'„ C&), 

A[a k ' + a(H(e x + H 2 AT) = A™ a + A™ [4.11a] 

b. for a transversely isotropic fiber ( C' rr = C (#) 

A^a = A™ a + A™ ~ [4.11b] 


3. Continuity of the radial stress a, at the fiber matrix interface, o'(a) = oT’(a) requires that: 


a. for a transversely orthotropic fiber 


<(Cjo + c'/.')a X ' - 1 + L/,: x + N'At 


<(C% + C) + A 2 m (C% - c”) --L - C™n™AT + C r % 

a 


or 


b. for a transversely isotropic Tiber 


[4.12a] 


A[(C^ + C f rr ) + C' rx r. x ~C f r] a'AT + 


= + c; 7) + A 2 m (C% - C) -i- + Cfr x - C”n”AT 


[4.12b] 
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4. Since there are no tractions applied at the outer boundary of the matrix o,(b) and t„(6) 


equal zero, the condition on a, requires that 


<(C™ + C) + /. 2 m -L (CZ - C) + C r % - C™a™AT - 0 [4.13] 

b 

The traction free condition on t„, requires that K : n Eqn. 4.8 equal zero. Thus, there are 
no shear stresses in the fiber or matrix. 


5. The Final condition for the case of pure thermal loading is that the net axial force P on the 
fiber-matrix composite be zero. This is expressed mathematically as 


p = 2*tj 0 6 o x «/r = 0 


[4.14] 


and for a transversely dhotropic fiber results in the condition 


<(C*'e + + <(0* 2 ~ ^ + y tC> 2 - 

a' +1 2 


+ L x a 2 >. x + -1 [A/ X a 2 -C™a™(b 2 - a 2 )]AF = 0 


For a transversely isotropic fiber (C,' r = C' M ), the form of Eqn. 4.14 is: 


[4.15a] 


y (C x ' rt + C xr )a 2 + <C X > 2 - a 2 ) 

- y [C> 2 -s 2 ) + C xx a 2 ]r x - -1 [cjuja 2 + c£a,V -a 2 )]AT - 0 


[4.15b] 


The constants 4{, Ap, 4". and t, are obtained by solving Eqns. 4 11-4.13 and 4.15 simul- 
taneously. It is noted that the axial loading case can be considered for a given axial strain 
c, or a given axial force P. In addition, the radial loading case may also be considered by ap- 
propriate modification of Eqn 4.13. 
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4.3 Results and Discussion 


Thermal stress distributions were determined for three types of fiber properties: (1) 
transversely isotropic (C„ = C w ); (2) circumferentially orthotropic (C„ < C M ); (3) radially 
orthotropic (C„ > C M ). The matrix was considered to be isotropic. The fiber and matrix prop- 
erties used for the calculations are given in Table 10. A uniform temperature increase of 
1°C was used for loading. Results for a fiber volume fraction V, = 0.623 are presented in Figs. 
93-95 and results for variable V, are presented in Figs. 96-98. 

Transversely Isotropic Fiber (C„ = C m ) 

The stress distributions for the case of a transversely isotropic fiber (Fig. 93) exhibit a 
uniform positive axial stress and uniform compressive hoop and radial stresses in the fiber. 
These distributions can be explained by examining the equation of the stress distribution in 
a transversely isotropic fiber (Eqn. 4.10a). Recalling that A\ = 0 in order to eliminate singular 
w displacements at r = 0 allows Eqn. 4.10a to be restated as 

of = A{(C' + Cl) + Clc x ~ cl.nl AT [4.16] 

It is evident from Eqn. 4.16 that the fiber stresses are independent of radial coordinate. The 
stress distributions in the matrix are described by Eqn 4.10c. A relatively large, positive hoop 
stress is present in the matrix. The hoop stress attains a maxim:'.",. at the fiber-matrix inter- 
face (Fig. 93). The decay in the magnitude of the radial and hoop matrix stresses with radial 
coordinate can be explained by examining the A" term in Eqn. 4 10c. The stresses decay as 
a function of 1/r*. The axial stress is constant in the matrix because C,» = C” when / = x Thus 
th» term in Eqn. 4.10c is zero and there is no axial stress dependence on the radial coor- 
dinate. For these strer distributions, fiber failure (if present) is expected to be an axial tensile 
fracture 
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Table 10. Fiber and matrix properties 


Radially Circumferentially Transversely Matrix 

Orthotropic Fiber Orthotropic Fiber Isotropic Fiber 


E, MSI(GPa) 

32(220) 

32(220) 

32(220) 

5(34.5) 

£„ MSI(GPa) 

4(27.5) 

32(220) 

4(27.5) 

5(34.5) 

E, MSI(GPa) 

32(220) 

4(27.5) 

4(27.5) 

5(34.5) 

Vrf 

.20 

.20 

.20 

.12 


.25 

.25 

.20 

12 


.025 

.25 

.25 

.12 

ct,(10 50 C) 

0.28 

0.28 

0.28 

1.11 

a«(10 *°C) 

5.56 

0.28 

5.56 

1.11 

a, (10 50 C) 

0.28 

5,56 

5.56 

1.11 
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Figure 93. Thermal stress distribution in a composite with a transversely orthotropic fiber. 
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Radially Orthotrooic Fiber (C„ > C*) 

Figure 94 shows the thermal stress distributions for the case of a fiber with radial 
orthotropy. All three components of stress are positive and singular at the center of the fiber 
These distributions can be explained by examining Eqn. 4.10b. For radial orthotropy 
C„ > C m and X{ is less than unity. Defining 5 = (1 - Xt). and recalling once again that 
A| = 0. Eqn. 4.10b can be written 

c' = 4(4 + cX)js- + Llz x + n!m [4.17] 

For the example problem considered here X{ = 0.343 and 5 = 0.657. Therefore, a stress 
singularity of order 8 exists at r = 0 for X( < 1.0. It should be noted that 5 is a function only 
of the fiber properties C' and CJ,. The term A((C£ + C£X$), which defines the strength of the 
singularity, is a function of fiber and matrix properties, fiber volume fraction, and boundary 
conditions. 

A similar singularity was shown by Lekhnitskii to exist in an anisotropic disk under radial 
compression. 4 * Lekhnitskii also notes that for the limiting case r = 0 there is no difference in 
material properties between the randG directions; therefore, the fiber must be transversely 
isotropic at r = 0. Such a condition precludes the existence of a mathematical singularity at 
the center of the fiber. (This point is also mentioned in a later paper by Olson and Bert.) 50 
However, the singular nature of the stress distributions as r approaches zero remains valid 
for the actual case. 

The potential failure mode of the fiber can be addressed by examining the relationship 
between the orientation of the graphite crystals in the fiber and the mechanical properties of 
a graphite crystal. For a radially orthotropic fiber, the basal planes of the graphite crystals 
Figure 90, which exhibit maximum strength are oriented parallel to the radial and axial di- 
rections “ Thus, the direction of minimum strength is in the hoop direction. Therefore, fiber 
splitting due to o» is a potential failure mode for a radially orthotropic fiber which exhibits 
singular hoop stresses. 
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Circumferentially Orthotropic Fiber (C M > C„) 

Stress distributions for circumferentially orthotropic fibers are shown in Figure 95. The 
stresses in the fiber are governed by the reduced form of Eqn. 4.10b. which is now written in 
the form 

of = /\((C4 + C^)r^ _1 + L* e x + N, f AT [4.13] 

and it is noted that >.{ - 1 > 0. Comparison of the stress distributions in Figure 94and 
Figure 95 (or comparison of equations 4.17 and 4.18) shows that the distributions in fibers with 
circumferential orthotropy (Figure 95) are completely different from those in radially 
orthotropic fibers (Figure 94). The distribution of axial and hoop stresses in circumferentially 
orthotropic fibers varies uniformly from compression along the centerline (r = 0) to tension 
at the fiber-matrix interface. The radial stress is compressive throughout the fiber. For this 
case A.; = 2.876; therefore, Eqn. 4.18 reveals that the stresses have a power function distrib- 
ution. The matrix exhibits compressive radial and axial stresses, but positive hoop stresses. 
All matrix stresses are relatively small in magnitude. 

In the circumferentially orthotropic fiber, the basal planes are oriented parallel to the 
axial and circumferential directions. Thus, the directions of maximum fiber strength are ex- 
pected to be in the axial and hoop directions with minimum fiber strength in the radial direc- 
tion. Therefore, the maximum tensile stresses are in the directions of maximum strength. 
The direction of minimum strength is under compression throughout the fiber. For this type 
of fiber orthotropy. failure, should it occur, would be expected to be via fiber splitting at the 
fiber-matrix interface. It is interesting to note that under a uniform temperature decrease the 
signs of the stresses will change, resulting in a positive radial stress. In this case the maxi- 
mum tensile stress will be in the direction of minimum strength. Therefore, the fiber may be 
more likely to fail during a cooling cycle than a heating cycle. 

It is noted that the solution presented here for a fiber in an isotropic matrix is quite dif- 
ferent than that of Chen & Diefendorf for a single fiber. 51 
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Figure 95. Thermal stress distribution in a composite with a circumferentially orthotropic fiber. 
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Influence of Fiber Volume Fraction 


The influence of fiber volume fraction on the distribution of thermal stresses is demon- 
strated in Figure 96 through Figure 98. Figure 96 shows the results for a transversely 
isotropic fiber. Figure 97 for a fiber with radial orthotropy, and Figure 98 a fiber with 
circumferential orthotropy. The fiber volume fraction was varied by changing the thickness 
of the matrix layer surrounding the fiber and holding the fiber radius constant, results are 
presented for fiber volume fractions in the range 0.391-1.0. 

These figures show that the axial component of stress is a function of fiber volume frac- 
tion for all three types of fiber microstructure. This is a direct consequence of the equilibrium 
requirement of zero axial force for pure thermal loading. Equilibrium must always be satisfied 
regardless of material properties. A somewhat surprising result is the fact that the distrib- 
utions of radial and hoop stresses are essentially independent of fiber volume fraction for both 
types of transversely orthotropic fibers considered (Figure 97b-c and Figure 98fc-c), but the 
distributions of these two stress components varies considerably with fiber volume fraction for 
the case of a transversely isotropic fiber. These results can be explained by considering the 
equations for the stress distributions in each type of fiber (Eqns. 4.16-4.18). 

In the transversely orthotropic fibers, A{ is relatively independent of volume fraction, dif- 
fering by less than 4 percent in the volume fraction range 0.391-1.0. In contrast, e, differs by 
more than 30 percent in the same volume fraction range. Detailed examination of the 
equations shows that the hoop and radial stresses are a strong function of the term containing 
Af and a weak function of the term containing e, since (C4 + A.(C') is large and L, is small for 
/ = rorO. Thus, it can be concluded that in transversely orthotropic fibers the hoop and radial 
stresses are relatively independent of fiber volume fraction. The axial stresses in the 
transversely orthotropic fibers, however, show a greater dependence on fiber volume fraction. 
In this case ( C C[ r ) is smaller and L, is more than two orders of magnitude greater than 
the corresponding terms mentioned above. Consequently, the axial stresses are a stronger 
function of fiber volume fraction through the term associated with axial strain. In transversely 
isotropic fibers, however, both Af and e, are strong functions of fiber volume fraction. A' and 
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e, vary by approximately 23 and 40 percent, respectively, in the fiber volume fraction range 
0.391 - 1.0. Thus, all stresses will be greatly influenced by fiber volume fraction. 

An equally surprising result is that the stresses in transversely orthotropic fibers are 
non-zero for a fiber volume fraction of 1.0, which corresponds to a fiber with no matrix sur- 
rounding it. The physical explanation is that as the fiber expands radially it also expands in 
the hoop direction; however, if there is a mismatch in the thermal expansion coefficients in the 
radial and hoop directions the expansion in the hoop direction can't compensate for the radial 
expansion. Thus, an internal constraint exists which gives rise to internal stresses. It should 
be noted that the presence of non-zero stresses in a single fiber is a function of the mismatch 
in radial and hoop thermal expansion coefficients only and not a function of material stiffness 
coefficients. 
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Figure 96. Effect of volume fraction on thermal stress distributions In a composite with a 
transversely isotropic fiber. 
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Figure 97. Effect of volume fraction on thermal stress distributions in a composite with a radially 
orthotrop'c fiber. 
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5.0 Conclusions 


A carbon-carbon composite manufactured from an 8-harness satin woven graphite cloth 
and a phenolic resin precursor was investigated both experimentally and numerically. The 
experimental portion of the study consisted of performing an out-of-plane tensile test in a 
scanning electron microscope and videotaping crack propagation. Additional scanning 
electron microscope studies were performed in order to determine the failure modes. The 
numerical portion of this study consisted of developing a finite element model and performing 
analyses on a two-layer woven fabric laminate for in-plane, out-of-plane, and thermal loading 
in order to characterize stress states in the laminate as a function of geometric parameters. 
Additionally, an elasticity solution was presented which investigates the influence of fiber 
anisotropy on thermal stress distributions in composites. A summary of the results of these 
studies are presented below. 

5. 1.1.1 Experimental 

1. Failure due to out-of-plane loading was interlaminar, occurring at the interface between 
adjacent plies. 
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2. Cracks propagated between the fiber-matrix interface. The cleanliness of fiber surfaces 
subsequent to failure indicated low fiber-matrix bond strength. 

5. 1.1. 2 Finite Element Analyses 

Out-of-Plane Loading 

.. The interlaminar normal stress n 2 is the dominant stress, and is relatively independent 
of stacking sequence, axial coordinate, undulation offset ratio, and undulation aspect ra- 
tio. 

2. The interfacial a, and components of stress are strong functions of stacking sequence 
and undulation aspect ratio. Additionally, o x and are strong functions of the undulation 
offset ratio R for R< 1, but relatively independent of R for R 2 : 1. The interfacial o y is 
relatively independent of geometric parameters. 

3. The best approach to increasing the out-of-plane tensile strength of a 2-D carbon-carbon 
laminate is to increase the fiber-matrix bond strength This is based on the following 
factors: (1) experimental observation of failure revealed that failure occurred between 
adjacent plies and that cracks propagated along the fiber-matrix interface; (2) finite ele- 
ment analyses indicates that the stress state between adjacent plies is relatively inde- 
pendent of material geometry. 

In-"1ane Loading 

1. Large interlaminar stresses are predicted in the region of undulation. This may explain 
the delamination observed in woven fabric graphite/polyimide composite by Wagnecz ”. 


Conclusions 


166 



2. A high in-plane normal stress a„ up to 170% of the far-field value, is predicted for the 
region of undu'ation, and is due to the bending of the tows. 

3. The presence of undulations reduces the in-plane modulus E, 6-18 percent from that of a 
comparable (0/90), laminate. 

« 

4. The in-plane modulus £ x is essentially independent of R for R £ 0.5. 

5. The predicted values of E , and E y in the (F/W ? /F) laminate differ by approximately 10 per- 
cent. 

6. All components of stress are affected by stacking sequence and undulation aspect ratio 
£1 The predicted values tend towards CLPT as the undulation aspect ratio increases. 

7. All components of stress are relatively independent of undulation offset ratio R for 1. 

Thermal Loading 

1. Large in-plane normal stresses are predicted throughout the laminate for thermal loading. 
The maximum predicted values of a, and c y are 25 psi/°F and 10 psi/°F, respectively. 

2. Interlaminar stresses (a, and t„.) on the order of 3-5 psi/°F are predicted for the region 
of undulation 

Finite Element With Variable Material Properties 

1. A finite element has been developed which allows the material properties to vary with 
element curvature. This allows curved domains to be modeled without mismatches in 
material properties at element boundaries and hence reduces the number of elements 
required to obtain an accurate solution. 
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5.1 .1.3 Effect of Fiber Anisotropy on Thermal Stresses 


1. The distribution of thermal stresses in a fiber reinforced composite material is affected 
significantly by the microstructure of the fiber. 

a. If the fiber exhibits radial ortholropy, the distributions of all three components of 
normal stress exhibit a singularity of type r~* . where the order of the singularity is 
a function of the radial and circumferential stiffness coefficients of the fiber. 

b. For circumferentially and transversely orthotropic fibers there is no singularity in the 
stresses. 

2. Fiber volume fraction has essentially no influence cn the radial and hoop stresses in 
orthotropic fibers. The axial stresses in orthotropic fibers and all three components of 
normal stress in transversely sotropic fibers are a function of fiber volume fraction. 

3. Singia fibers exhibit non-zero stresses when there is a mismatch in the radial and 
circumferential thermal expansion coefficients. 
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Appendix A. Contour and Interfacial Stress 
Distributions 


The following figures show the stress contours for the (F/W^F),, R — 1 and 4 laminates for 
in-plane loading e„ = 0.1% and out-of-plane loading e r = 0.1%. Additionally, interfacial stress 
distributions are presented for those (F/Wj/F),, (F/Wj/F)„. and [F/W) z laminates that were not 
discussed in detail Chapter 3. 
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Figure 101. Contour plot o 
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Figure 104. Contour plot of in (F/Wj/F^ laminate. R=4, e z 
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Figure 105. Contour plot of O j in (F/W ? /F) ( laminate. R=4, =0.1%. 
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a) Fill Tows 



A = 4psi 

b) Warp Tows 


Figure 106. Contour plot of in (F/W ? /F) ( laminate. R=4, =0.1%. 
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b) Warp Tows 


t-'igure 109. Contour plot of in (F/W } /F) ( laminate. R = 1, r.^ = 0.1%. 
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Figure 112. Contour plot of In V 'F/W J /F) ( laminate. R = 4, = 0.1%. 
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Figure 114. Contour plot of t >; In (F/W 2 /F) ( laminate. R =4, e* = 0.1%. 
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Figure 118. Interfacial C t and t |t in (F/W 2 /F) ( laminate. R=4, = 0.1%. 
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Figure 125. Interfacial O i and in (F/Wj/F^ laminate. R = 1. r.^ = 0.1%. 
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Figure 127. Interfacial o * and in (F/W ? /F) ( laminate. R = 4, r.^ — 0.1%. 
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Figure 129. Interfacial 0 and 0 in (F/VV „/F) laminate. R = 1,r. =0.1%, 
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Figure 130. Interfacial a and t in (F/W /F) laminate. R = 1,e =0.1%. 
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Figure 131. Interfacial CT and o in (F/W /F) laminate. R = 4, r. = 0.1%. 
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Figure 132. Interfacial o and x in (F/W/F) laminate. R=4, F. =0.1%. 
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Figure 133. Interfacial and in (F/W) ? laminate. R = 0.5, r. < = 0.1%. 
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Appendix B. influence of Boundary Conditions on 
Finite Element Analysis 


One characteristic of a two-layer woven fabric composite is that the geometry is periodic 
in x (Figure 17). Thus it can be divided into repeating unit cells (Figure 20a). Ideally the re- 
sponse of an entire laminate can be modeled using a single unit cell. This simplifies the 
problem and reduces effort and computational expense. A disadvantage of modeling a single 
unit cell is that the response of the laminate at the boundaries of the unit cell (A-A and B-B 
in Figure 20a) is unknown, yet any analytical or numerical analysis requires the imposition 
of some form of boundary condition at those locations. 

In this study, a unit cell of a woven fabric composite was modeled using the boundary 
conditions illustrated in Figure 20. For both out-of-plane and thermal loading the ends of the 
unit cell (x=0, x = L) were constrained to remain straight and vertical. For in-plane loading 
the left end of the unit cell (x=0) was constrained to remain straight and vertical and the right 
end (x = L) was given a displacement corresponding to an axial strain e, = 0.1%. In this ap- 
pendix a comparison of the interfacial stress distributions for one and two unit cell models are 
made in order to test the validity of the boundary conditions stated above The comparison 
is made for both in-plane and out-of-plane loading. No comparison is made for thermal 
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loading because the boundary conditions on the ends of the unit cell are the same for both the 
thermal loading and out-of-plane loading cases. 

The analysis was done by building one unit cell and two unit cell finite element models 
of *hs (F/Wj/f)^ R=4 laminate. The one unit cell model contained 452 elements and the two 
unit cell model contained 912 elements. This particular laminate was chosen because the 
undulations, which exhibit a complex stress distribution, are close to the ends of the models. 
Thus, this laminate provides a worst case test of the effect of the boundary conditions of the 
stress distribution in the laminate. The boundary and loading conditions for out-of-plane and 
in-plane loading corresponded to those illustrated in Figure 20b and c, respectively. 


B.1 Out-Of-Plane Loading 


The interfacial stress distributions of the one and two unit cell models subjected to out- 
of-plane loading are compared in Figure 139 through Figure 142. The in-plane stresses o, 
and o y show no variation in stress between the two models, nor does the interlaminar stress 
o, . The distribution (Figure 142), however, does indicate a difference between the one and 
two unit cell models at region B. At £ = 14. is zero in the one unit cell model but is ap- 
proximately 7 psi in the two unit cell model. This difference is a direct result of the constraint 
on the one unit cell model that the laminate remain straight and vertical at x = L (£ = 14). 
This constraint prevents shear deformation of the laminate at the boundary of the laminate 
(x = L) and thus requires t„ = 0. The two unit cell model is not constrained at £ = 14 and is 
thus allowed to shear. 

The effects of the boundary conditions are also evident at regions A and D. At £ = 0 
and £ = 30 . t„ is zero. From the distribution exhibited in region B, it can be concluded 
that should more than two unit cells of the laminate be modeled the distribution in region D 
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would be the same as that of the two unit cell distribution at B and that the true r„ distribution 
at region A should be the same as that at region C. 

Away from the unit cell boundaries there is no difference in the t„ distributions between 
the one and two unit cell models, nor is there a difference in the maximum or minimum value 
of in the laminate. 


B.2 In-Plane Loading 


The interfacial stress distributions of the one and two unit cell models are compared in 
Figure 143 through Figure 146. All components of stress show no variation with the number 
of unit cells modt i*. This i attributed to two factors. First of all, under in-plane loading the 
right side of the model (x = L) is under load. Therefore, no artificial constraint exists at that 
location Secondly, under in-plane loading the interfacial stress distributions rapidly converge 
toward CLPT away from the undulation, ai:d thus the boundary condition at x=0 has no effect. 


B.3 Conclusions 


The results discussed above indicate that the boundary conditions on the one unit cell 
model have no influence on the stress distribution for the in-plane loading -nse. For the out- 
of-plane loading case only the interlaminar shear stress is affected, and the effect is rela- 
tively small and limited to a small region near the boundary. Tnus, it can be concluded that 
the woven fabric laminates analyzed in this study can be accurately modeled with a single unit 
cell. 
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Figure 139. Comparison of interfacial cf^ distributions for one and two unit cell FEM models under 
out-of-plane loading: (F/W ? /F) ( laminate, R— 4, e^-O.1%. 
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Figure 140. Comparison of interfacial o distributions for one and two unit cell FEM models under 
out-of-plane loading: (F/wVF) ( laminate, R-4, r.^-0,1%. 
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Figure 142. Comparison of interfacial t distributions for one and two unit cell FEM models under 
out-of-plane loading: (F/W*/F) ( laminate. R-4, r.^-0.1%. 
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Figure 144. Comparison of interfacial a distributions for one and two unit cell FEM models under 
axial loading: (F/W ? /F) ( laminate, R-4, r^-0.1%. 
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Figure 146. Comparison of interfacial r distributions for one and two unit cell FEM models under 
ixlal loading: I F W^/F^ lo.ninate, R - 4, e > - 0. 1 %. 
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Appendix C. NVT3D User's Guide 


C.1 Introduction 


NVT3D is a linear elastic finite element code developed under the assumption of gener- 
alized plane strain. The assumption of generalized plane strain allows the user to study the 
response of an infinitely long body of arbitrary cross-section to various thermal-mechanical 
or hygro-mecharical loadings with the restriction that all quantities except axial displacement 
are independent of the axial coordinate z. This appendix is the user's guide to NVT3D. It 
contains the theoretical development of the finite element model for generalized plane strain 
and presents an overview of the program's capabilities, the input requirements, and output 
capabilities. 

The formulation and assembly of the element stiffness matrices is based upon the method 
of Reddy 11 The method of storing and solving the resulti g system of equations is based 
upon the finite element program STAP (STatic Analysis Program) developed by Bathe and 
Wilson **. 
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C.2 Generalized Plane Strain Finite Element Formulation 


The most general displacement field for generalized plane strain problems can be written 

u - mxy) [C.la] 

v = V[xy) [C.lb] 

w 0 = HT(x,y) + z z z [C.lc] 

The coordinate system used by NYT30 is shown in Figure 147. The x.y and z axes represent 
the global coordinate system while the material principal directions are represented by the 
..2 and 3 axes. The region to be modeled oy finite elements is then a typical cross-section, 
which occupies the x-y plane. 

For generalized plane strain the reduced from of the equilibrium equations are 


CO. 

+ 


rx 

*V 

‘ ^xy ^°y 

cx 

<5y 

4 

ci yz 


cx cy 


= 0 


= 0 


and the reduced form of the strain-displacement r ations are 


[C.2a] 


[C.2b] 


[C2c] 


[C.3a] 
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Figure 147. Definition of Material (1-2-3) and Global (x-y- 2 ) Coordinate Systems Used in Finite 
Element Formulation 
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[C.36] 


Ky dy 


„ _ dw 

v IF 


= — + ~ 
ay cx 


The linear stress-strain relationships are given by: 


ty ~ OyAT 


a xy ~ a xy& T \ 


[C ] is the material stiffness matrix in terms of the global (x-y-z) coordinate system given by: 


[c] = C tJ-’CcM 


[C] is the material stiffness matrix in the material (1-2-3) coord’nate system and is given by: 
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where: 


[C] = 


'12 "13 


C22 Cj 3 


G 11 ~ 


1 - Vo-aV-y 


-11 "12 


c„ — 


_ v 12 + v 13 v 32 , 


^ _ v 13 + v 12 v 23 , 

C 13 ^ 1 


C 22 - 


1 - v„v. 


C23 ~ 


_ v 23 + v 21 v 13 , 


1 - V 12 V 21 


G 44 “ G 23 


C 5$ “ G 13 


G 66 ~ G 12 


A = 1 - v 12 v 21 - v 13 v 31 - - v 12V23 v 31 - v 21 v 13 V 32 


and the transformation matrices C^] and [T z ] for a rotation 0 about the 2 axis are given by: 


m 2 n 2 0 0 


n 2 m 2 0 0 


-2 mn 


W- 


0 1 


0 m 


0 n 


-mn mn 0 0 


0 m 2 -n 2 
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Cr 2 ] = 


2 2 
m n 


n 

0 

0 

0 


m 

0 

0 

0 


-2 mn 2 mn 


0 

0 

1 

0 

0 

n 


0 

0 

0 

m 

n 

0 


0 

0 

0 

—n 

m 

0 


mn 

-mn 

T 

0 


2 2 
m - n 




where m = cos 0 n = sir. 0 

Substituting [C.3a-e] and (C.4J into [C.2a-cl yields the governing equations for the defor- 
mation of a plane elastic body under generalized plane strain: 

■Jr[Mlr " aAT ) + Mi “ QyAr ) + ^ ” aA1) 


+ nt + & - °^ r )] + iNt - °^ t ) + m - °’ ar ) tc8a] 

+ ^s(«* - 0^7) + + jL- a xy ATjj = 0 

■J^Mi “ a * AT ) + Mi ~ Q / AT ) + ^ ~ 


+ Mi + f - a -M] + i[Mi - a ^ T ) + Mi - [c 8b] 


C»(c x - a^T) + C 26 (^ + -g- ~ a xy AT)j = 0 
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iMrhiKf]' 0 


[C.8c] 


However, an additional equation is needed to govern the axial strain (eJ response, which is 
a variable in the generalized plane strain formulation. Thus, the net axial force F t must be 
specified. This is expressed mathematically as 

fojdxdy = F z [C.9a] 

In terms of the displacements it is of the form: 

f{Mlr ' “«") + Mt “ a ^ r ) 

[C.9b] 

+ C33(e z - *£T) + + jL - a xy Ar)^dxdy = F z 

Taking the variational form of [C.8a-c] and (C.Stl, and separating them into terms associated 
with toe primary variables u. v, w, and £ r yields the following finite element stiffness matrix: 


[«’>] * * C,jz- + C,A — + c3- i? 'L 

Jy 1 dx dx 16 dx dy 16 by dx 66 by dty 

o' 3 ] - f (c,.^L + c,J*- ^ i ? I + 

‘ 5x dy 16 dx dx 2er dy dy 66 dy dx 


' 16 17 

dy 

r <V/ 

dy t 

16 dx 

dx 

r 

dy. 

66 dx 

dx 

r 

•'W _ ST 

dvj 


[Kl4] = j(lT Cl3 + ‘^ LC3e ) dXdy 


[C.lOa] 

[C.lOb] 

[C.IOc] 

[C.IOd] 

[C.lOa] 
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ck24] ~ {( it 0 * 5 + lfa) dxdy 


[C.lOf] 


[K 34 ] = 0 

tC.IOg] 

[K 44 ] = C33 

[C.lOh] 

[<] = [k; 2 ] [K 13 ] = [K 23 ] = [K 32 ] = [K 31 ] = 0 

[C.IOi] 


The equations on the secondary variables are 


{f1} = f{^ c,iaxAr + c '* a ** T + c '* aAT + c '&*y* T ) 

[C.lla] 

+ '^ L ( C 1 6<1 x Ar + C 26 a y A7 ‘ + p36 a z A7 ’ + C 66 a *jA r )\ A7 ri /l + J r V,f//T 


{p2} = j{^(^6«xAr + CtfiyAT + + CvPtyAT) 

[C.llb] 

+ ^(^ a * AT + ^2°yAT + C^ az AT + C 26 a xy AT)\ArdA + ^t/JT 


{F 3 } = J rV( ^r 

[C.llc] 

II 

[C.lld] 

and represent the nodal forces and equivalent thermal loads. 
The resulting system of linear equations is of the form: 
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[K 11 ]^} + [K 12 ]{v> + [* 13 ]{w} + [K 14 ]{e 2 > = {F 1 } 


[K 21 ]{o> + [K 22 ]M + [K 23 ]{w> + [K 24 ]{e z } = {F 2 } [C. 12] 

[K 31 ]{u> + [K 32 ]^} + [K 33 ]^} + [K 34 ]^} = {F 3 } 

[K 4, ]{u> + [K 42 ]M + [K 43 ]{w> + [K M ]{e z } = {F 4 } 

It should be noted that e* is not a nodal degree of freedom, but is constant over the entire 
model. Therefore, e z as a variable results in the addition of only one equation to the global 
stili.iess matrix. 


C.3 Element Geometries 


Three elements are available in NVT3D. They are the 6-node triangular isoparametric 
element, the 8-node quadrilateral isoparametric element, and the SP-8 (SPecial 8-node) ele- 
ment. The geometries of the elements along with the node numbering and Gauss point 
numbering schemes are illustrated in Figure 148. 

The 6-node triangular and 8-node quadrilateral elements are standard orthotropic ele- 
ments. The SP-8 element n&s the same geometry and node numbering as the standard 
8-node quadrilateral element, but in the SP-8 element the material stiffness matrix [c] is 
evaluated at each Gauss point in the element. This feature allows the material properties to 
vary within an element It is useful where there is significant curvature in a model (e.g. 
axisymmetric geometries). This element allows curved geometries to be modalled more ac- 
curately with fewer elements. 
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Figure 148. Available Element Geometries 
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C.3.1 Formulation and Use of the SP-8 Element 


In practice, if a finite element is of a general two-dimensional shape, the integration of the 
IC.11] must be done numerically. The numerical integration is done by mapping the element 
into a ’master* element which is a square for an 8-node element (Figure 149a). This method 
is convienent because the master element has constant limits of integration. The mapping 
functions are: 


X = 2x ( v, 

[C.13a] 

y * Zy.vi 

[C.13b] 


where x, and y, are the nodal coordinates in the global (x-y) system and y, = y,(!;, q) are the 
finite element interpolation functions. 

Taking the derivative of (C.13] with respect to 5 yields the following equations: 




Dividing [C.14b] by [C.14a] yields the form of the equation to determine the slope of any point 
in the global coordinate system with respect to the \ axis: 
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If the element material properties are described in the S-tj coordinate system then the trans- 
formation angle 0 between the material coordinate system and any point in the general ele- 
ment can be determined: 


0 = tan 


’(f) 


[C.16] 


Hence, the material stiffness matrix [C] can be calculated at any point in the element using 
[C.6-7], 

This method is much simpler, though less general, than the parametric cubic modeling 
approach used by other researchers, which requires modification of the interpolation functions 

58 59 60 


C.3.2 Implementation of SP-8 Element 


To use the SP-8 element correctly, the correct element type must be specified (see INPUT 
REQUIREMENTS) and the material properties must be input such that the 1 and 2 directions 
in the material coordinate system correspond to the % and q directions, respectively, in the 
element. The best way to do that is to recall that the node numbering scheme is such that the 
line (or curve) formed by nodes 1, 5, and 2, lies along the \ direction. As such, the following 
relationships hold. 


£ i = E (, £ 2 = £ n £ 3 = E z 
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a) 8-Node Element of General Shape Mapped Into Master Element 



b) Transformation of Material Properties With Change In Element Shape 


Figure 149. Mapping and Coordinate Transformation of SP-8 Element. 
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v 23 * v V V ‘3 “ V V « = 

G 23 “ G m G 13 = G 12 = G 4m 

°1 ~ °i ! a “> = ^ c 3 ~ a z 


Therefore, the user inputs the elastic constants in the t-q coordinate system and NVT3D cal- 
culates the transformation anyle 6 anu subsequently [C] at each Gauss point in the element. 


C.3.3 Verification of SP-8 Element 


The ai..uracy o f the SP-8 element was verified by comparing the finite element solution 
for a cylindrical^ orthctropic fiber under unifor~i thermal load with the elasticity solution 
presented in Chapter 4. ne results were also compared with a finite element model using 
standard 8-node elements. 

Two nme-elrment. quarter-symmetry finite models of a circumferentially o^hotropic fiber 
under a uniform thermal toart AT — ? C C were analyzed (Figure 150). Model A contained three 
6-node triangular elements and six standard 3-node elements The average transformation 
angle 0 of i arh element was determined and entered into the input file. Model B contained 
three 3-node elements and six SP-8 elements The average transformation angie h was en- 
tered for the 6-node triangular element only. For the 3P-8 elements the materia' properties 
were input suta that tha elastic constants in the 1-2 system corresponded io those of the c-q 
coord mite c ystem 

T t-e hoop radial, and axial stresses (c„, x, and cr s ) of the tnree soiut'or.s are compare.' 
ic r eac. Gauss point in the mode! >n Figure 151. Figure 15? and Figure 153. respec vely 
Smcr ‘he problem -s rnmetr e tn^ stress distribn’ - n at all ihe ~ points ; n the mesh 
houlo be equal The resu : ts indicate excellent agreement betv .. . the eicticity solution and 
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the finite element solution using the SP-8 element. The difference between the solutions is 
less than 5% for all three components of stress. The finite element solution using the stand- 
ard 8-node elements differs from the elasticity solution by up to 13% for o, and by up to 35% 
for the o v and a, components of stress. 
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Figure ISO. Finite element models of circumferentially orthotropic fiber under uniform thermal 
load. AT = 1*C. 
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Figure 152. Comparison of FEM and exact solutions to O f distribution for a fiber under thermal 
load. 
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C.4 NVT3D Input Requirements 


NVT3D was developed assuming that the input file would be submitted from a CRT and 
not typed on cards. Thus, the input structure has been designed such that there is room cn 
each card for description(s) of the variable(s) that are read from that card. In this section a 
description of each type of card image in the input file is given. The required format 
card image is given in Section C.4. 

• Title (Cards 1 and 2) 

■ The title consists of two 80 character lines which will be printed on the output flies 
and on the plots. 

• Number of Element Groups (Card 3} 

« On this card the user specifies the number of element group? in the mesh. An ele- 
ment group consists of a set of elements, ea^h having the same element Type, the 
same material properties, and the same orientation 0 (see Fig. 1). 

• Number of Noues in the Mesh (Card 4) 

• Nodal Coordinate X and Y Scale Factors (Card 5 & 6) 

■ To ease input the user can specify scale factors which independently scale the 
magnitude of the x and y nodal coordinate values. 

• Data Check Flag (Card 7) 

• This option allows the user to check the input without proceeding to the solution. 
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• Print Nodal Coordinates (Card 8) 


• This option allows the user to inhibit the output of the nodal coordinates. This opt'on 
is helpful in reducing the amount of output if multiple jobs are run using the same 
mesh. It should be noted that the nodal coordinates are always printed if the job is 
in the data check mode. 

• Plotter Flags (Cards 9 & 10) 

■ On this card the user sets the flags for plotting the undeiormed and deformed 
meshes. 

• Element Range to Plot (Card 11) 

• Sometimes it is desirable to plot just a portion of the finite element mesh. This option 
allows the user to specify the range of elements to be plotted. 

• X and Y Offsets for Mesh Plot (Cards 12 & 13) 

■ This option allows the user to shift the X and Y coordinates of the plots. It may be 
useful when using the 'Element Range to Plot' command (Card 11). 

• Mesh Plot X and Y Scale Factors (Cards 14 & 15) 

■ This option allows the user to scale the size of the plot. This option is especially 
useful when used with the Range to Plot card (Card 11) to plot very small elements. 
The y-coordinate scale factor (YMSCAL) must be set such that the maximum y- 
coordinate in the mesh will be scaled to a value less than or equal to nine inches. 

• Maximum U and V Displacements (Cards 16 t\ 17) 
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This option allows the user to specify in inches the maximum u and v displacements. 
All of the smaller displacements are scaled proportionally to this value. 


• Plane Strain/Plane Stress Flag (Card 18) 

• NVT2D: Specifies whether analysis is plane stress or plane strain. 

■ NVT3D: Not used. 

• Plate Thickness (Card 19) 

- Specifies the thickness of the plate being analyzed. 

• Thermal Analysis Flag (Card 20) 

■ This flag determines whether thermal (or hygro) analysis is implemented. 

• Temperature Change (Card 21) 

■ This card contains the temperature for thermal-mechanical loading or the moisture 
content for hygro-mechanical loading. 

• Number of Specified Nodal Degrees of Freedom (Card 22) 

■ The user inputs the nunber of specified displacements in the problem. 

• Number of Specified Boundary Forces (Card 23) 

* The user inputs the number of specified forces in the problem. 

• Number of Constrained Nodes (Card 24) 
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• The user has the option of specifying that certain nodes have the same, though un- 
known, displacements (e.g. an edge of r. model uv.'st deform such that it remains 
straight). Examples of the use of this feature and the finite element formulation used 
to implement this feature can be found in Branca * Crane “and Adams 57 . This card 
specifies the number of nodes that are constrained to have the same displacement. 

• Axial Load Flag (Card 25) 

■ Specifies whether axial strain e z or average axial force F x is applied. 

• Magnitude of Axial Load (Card 26) 

• The user inputs the magnitude of the applied axial strain z t or the applied axial force 
F t . 

• Gauss Point Stress, Strain, ?nu Strain Energy Dumps (Cards 27-29) 

• These options allow the user to dump all of the stresses, strains, and strain energies 
at all of the Gauss points to separate output files (see Output Capabilities for more 
information). 

• Nodal Point Displacement Dump (Car 30) 

■ This option allows the user to d'imp the displacements of all of nodes to a separate 
output file. 

• Nodal Point Stress, Strain, and Strain Energy Dumps (Cards 31-33) 

■ These options allow the user to dump ail of the stresses, strains, and strain energies 
at all of the nodes to separate output files (see Output Capabilities for more infor- 
mation). 
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• Element Group Data (Cards 34-40) 


• Element Group Control Cards (Cards 34-3 9) 

A Card 34 is blank. 

a Card 35 contains the element type, the number of elements in the element group, 
the orientation angle of the elements with respect to the x-axis, and the param- 
eter delta (delta is. unused at tt ~ time). 

a Card 36 contains the elastic moduli E„ E z , and £ 3 . 

a Card 37 contains the Poisson's ratios v a , v 13 , and v, 2 . 

a Card 38 contains the shear moduli G 23 , G 13 and G 12 . 

a Card 39 contains the thermal expansion coefficients a„ a 2 , and a 3 or the moisture 
expansion coefficients p t , p 2 , and P 3 . 

• Element Connectivity and Print Card (Card 40) 

a This card contains the element connectivity matrix and output flags. This card 
is repeated for each element in the element group. See Fig. 2 for element 
numbering schemes. 

■ Cards 34-40 are repeated NEG times, where NEG is the number of element groups in 
the mesh. 

• Nodal Data Information (Card 41) 
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• The node number and node coordinates are specified on this card. Card 41 is re- 
peated NNOD times, where NNOO is the number of nodes in the mesh. 

• Specified Nodal Displacement Information (Card 42) 

• On this card the user identifies the node number, direction, and maQnitude of the 
displacement of each specified degree of freedom. Card 42 is repeated NSDF times, 
where NSDF is the number of specified degrees of freedom. 

• Specified Nodal Force Information (Card 43) 

• On this card the user specifies the node number, direction, and magnitude of each 
of the specified nodal forces. Card 43 is repeated NSBF times, where NSBF is the 
number of specified boundary forces. This card is omitted if NSBF equals zero. 

• Specified Nodal Constraint Information (Card 44) 

■ On this card the user specifies the node number and direction for each of the con- 
strained degrees of treedom. This card is repeated NCON times, where NCON is the 
number of constrained displacements. 


C.5 NVT3D Outnnt Capabilities 


NVT3D solves foi the noda displacements, stresses, strains, and strain energy densities 
for each element in the mesh. These values are calculated in . ie global (x, y, z) and principal 
material (1, 2, >' coordinate systems. 
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The user has the ability to tailor the output file. By setting the output flags on Card 41. 
he can specify the output for each Gauss point and nodal point in the element. The output 
points for the eight-node and six-node elements are shown in Fig. 3. 

To specify the Gauss point output for the eight-node element (six-node element) the 
Gauss point numbers for which output is desired are placed in columns 47-55 (47-53). For 
example, suppose output is desired for Gauss points 1. 3, 5, and 9. Columns 47-55 would 
therefore contain 135900000 . The order of the input does not matter. The input could be any 
of the following: 

195300000 

001090350 

000013590 

etc. 

If output is desired for all of the Gauss points in the element columns 47-55 would read 
123456789. If no output is desired for the element columns 47-55 would contain all zeros. 

To specify the nodal output for the eight-node (six-node) element the local node numbers 
for which output is uesi.ed are placed in columns 57-64 (57-62). The input format is the same 
as that of the Gauss points. 

It should be noted that the Gauss point output flags and the nodal point output flags are 
eacn read as one integer variable. This reduces execution time and the amount of memory 
required As a result, when the output flags are echoed to the output file a single zero will 
be printed in column 55 (53) if no Gauss point output is requested for the element. Similarly, 
if no nodal output is requested a single zero will be printed in column 64 (62). 

Frequently it is desired to have all of the displacements, stresses, strains, or strain en- 
ergy densities output in a format that is more conducive to post-processing. The output 
'dump' flags (Ca> "’7-33) allow the user to do this. Activating the appropriate Hags will cause 
separate output files to be generated that contain the displacements, stresses, strains, or 
strain energy densities. Each line of the output file will contain the element number, the output 
point nuinbei (node number or Gauss point number), the coordinates r ' the output point, and 
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the values of the displacements, stresses, strains, or strain energy densities (depending on 
which flag is activated). 


C.6 NVT3D Input Card Sequence 


Card 

1 

2 

3 

4 

5 

6 
7 


8 


9 


10 


11 


Format 

(20A4) 

(20A4) 

(T57.I2) 

(T57.I4) 

(T57.E10.3) 

(T57.E10.3) 

(T57I1) 


(T57.11) 


(T57.I1) 


(T57.I1) 


(T57.I3.1X.I3) 


Variable 

ITIT1 

ITIT2 

NEG 

NNOD 

XSCAL 

YSCAL 

NCHK 

0 = Data Check Mode 

1 = Solution Mode 


Contents 

Title 

Titie 

Number of Element Groups 
Number of Nodes in Mesh 
X-Coordinate Scale Factor 
Y-Coordinate Scale Factor 
Data Check Option 


NCPRT Nodal Coordinate Print Option 

0 = Un-activated (unless NCHK = 0) 

1 = Prints Nodal Coordinates 


MESHU Undeformed Mesh Plot Option 

0 = No Plot 

1 = Undeformed Mesh Plotted 

2 = Undeformed Mesh Plotted (Node & Element No.) 


MESHD Deformed Mesh Plot Option 

0 = No Plot 

1 = Deformed Mesh Plotted If NCHK = 0 


LOWPLT First Element Number Plotted 

IHIPLT Last Element Number Plotted 

= 399 Then All Elements > = LOWPLT Plotted 
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12 (T57.F84) 

13 (T57.F8.4) 

14 (T57.E10.3) 

15 (T57.E10.3) 

16 (T57.E10.3) 

17 (T57.E10.3) 

18 (T57.I1) 


19 (T57.E10.3) 

20 (T57.il) 


21 (T57.E10.3) 

22 (T57.I2) 

23 (T57.I2) 

24 (T57.I2) 

25 (T57.I2) 


26 (T57.E12.5) 


27 (T57.I1) 


28 (T57.I1) 


XOFF 

YOFF 

XMSCAL 

YMSCAL 

UMAX 

VMAX 

NPS 

0 = Plane Strain (NVT2D) 

1 = Plane Stress (NVT2D) 
Not Used By NVT3D 


X-Axis Offset of Plot 
Y-Axis Offset of Plot 
X-Coord. Scale Factor of Plot 
Y-Coord. Scale factor of Plot 
Max. U Displ. of Deformed Mesh 
Max. V Displ. of Deformed Mesh 
Plane Stress/Plane Strain Option 


THICK Plate Thickness 

NHEAT Thermal/Hygro Analysis Option 

0 = No Thermal/Hygro Analysis 

1 = Thermal/Hygro Analysis 


TEMP 

NSDF 

NSBF 

NCON 


Temoerature/Moisture Content 
Number of Specified Displacements 
Number of Specified Forces 
Number of Constrained Displacements 


NEZ Axial Load Flag 

0 = Axial Strain is Specified 

1 = Axial Force is Specified 

VEZ Magnitude of Axial Load 

= r. t if NEZ = 0 
= F t if NEZ = 1 

NGP(1) Gauss Point Stress Dump Option 

0 = Unar’ivated 

1 = Global Coordinate System Only 

2 = Local Coordinate System Only 

3 = Global and Local Coordinate Systems 


NGP(2) 


Gauss Point Strain Dump Option 
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0 = Unactivated 

1 = Global Coordinate System Only 

2 = Local Coordinate System Only 

3 = Global and Local Coordinate Systems 


29 (T57.I1) NGP{3) 

0 — Unactivated 

1 = Activated 

30 (T57.I1) NPD(1) 

0 = Una Pir ated 

1 = Activated 


Gauss R. Strain Energy Dump Opt. 


Nodal Displacement Dump Option 


31 (T57.I1) NPD{2) Nodal Stress Dump Option 

0 = Unactivated 

1 = Global Coordinate System Only 

2 = Local Coordinate System Only 

3 = Global and Local Coordinate Systems 


32 (T57.I1) NPD{3) Nodal Strain Dump Option 

0 = Unactivated 

1 = Global Coordinate System Only 

2 = Local Coordinate System Only 

3 = Global and Local Coordinate Systems 


33 (T57.I1) NPD{4) Nodal Strain Energy Dump Option 

0 = Unactivated 

1 = Activated 


34 Blank Card 

35 (T1S,I1.T36.I4.T48.F6.2.T62,F6.2) 

ITYPE Element Type 

1 = 6-Node Triangular Element 

2 = Not Used 

3 = 8-Node Quad. Element 

4 = SP-8 Element 

NUMEL Number of Elements in Group 
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ANGLE 

Angle of Elements With Respect 
To Global Coordinate System 


DELTA 

Not Used At This Time 

36 

(3(8X.E10.3)) PROP(1-3) 

Elastic Moduli £,, £,, and £ 3 

37 

(3(8X.E10.3)) PROP{4-6) 

Poisson's Ratios v„, v J3 , v, 3 

38 

(3(8X,E10.3)) PROP(7-9) 

Shear Moduli G a , G, 3 , and G I2 

39 

(3(8X,E10.3)) PROP(10-12) 

a,, a*. a 3 or P„ P*. p 3 

40 

(T2,6(1X,I4),T47,I7,T57,I6) 

Six-Node Triangular Element 


NOD(1-6) 

Element Connectivity Matrix 


KEYPT(I) 

Gauss Point Print Flags 


KEYPT(2) 

Nodal Print Flags 



or 


(T2,8(1X,I4),T47.I9.T57 I I3) 

Eight-Node Quadrilateral Element 


NOD(1-8) 

Element Connectivity Matrix 


KEYPT(I) 

Gauss Point Print Flays 


KEYPT{2) 

Nodal Print Flags 


Repeat Card 40 for each element in the element group. 
Repeat Cards 34-40 for each element group in the mesh. 

41 

(Free-Format) 

Input Nodal Coordinates 


IXY 

Node Number 


•'(IXY) 

X-Coordmate of Node IXY 


Y(IXY) 

Y-Coordinate of Node IXY 


Repeat Card 41 NNOD Times 


42 

(Free-Format) 

Input Specified Displacements 


ND 

Node Number 


NUDIR(ND) 

Direction of Specified Disp. 
(1 = u, 2 = v, 3 = w) 


VSD(ND) 

Magnitude of Specified Disp. 


Repeat Card 42 NSDF Times 


43 

(Free-Format) 

Input Specified Forces 


NF 

Node Number 
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NFDIR(NF) Direction of Specified Force 

(1 = u. 2 = v) 

VSF(NF) Magnitude of Specified Force 

Repeat Card 43 NSBF Times 
Omit Card 43 IF NSBF = 0 


44 


(Free-Format) 


NCD 

NCDIR(NCD) 


Input Constrained Displacements 
Node Number 

Direction of Constrained Displacement 
(1 = u. 2 = v. 3 = w) 


Repeat Card 44 NCON Times 
Omit Card 44 IF NCON = 0 


C.7 Implementation of NVT3D 


NVT3D was developed to run on and IBM 3090 operating under MVS or under VMBATCH. 


C.7.1 Implementation under MVS 

NVT3D may be executed under MVS by submitting the following files: 
NVT3D JCL1 
NVT3D FORTRAN 
NVT3D JCL2 
filename DATA 
NVT3D JCL3 
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C.7.1.1 NVT3D JCL1 


//Annnxxx JCB acct#.USERID,REGION = 1500K.TIME = 1 

/•PRIORITY IDLE 

/•JOBPARM LINES = 50.CARDS = 0 

//STEP1 EXEC FORI VCGV, PARM . FORT = 'NOSOURCE.NOSRCFLG' 
//FORT.SYSIN DD * 


where nnn is the user's output box number and xxx is the user's initials. The REGION and 
TIME statements may have to be changed depending on the size of the problem and the time 
needed for execution. 


C.7.1.2 NVT3DJCL2 


r 

//GO.FT06F001 DD SYSOUT = (A„1) 

//GO.PLOTLOG DD SYSOUT = A 

//GO.VECTR1 DD DSN = &&VECTR1 .UNIT = 3380, SPACE 

// DISP = v.f’ASS) 

//GO VECTR2 OD DSN =&&VECTR2, UNIT = 3330.SPACE 
// DISP = (.PASS),DCB = (BLKSIZE = 32002) 

//GO PLOTPARM DD DUMMY 
//GO.FT08F001 DD SYSOUT = (A.,2) 

//GO.FT09F001 DD SYSOUT = (A„3) 

//GO.FT11F001 DD SYSOUT = (A,.4) 

//GO.FT12F001 DD SYSOUT = (A,.5) 

//GO FT15F001 DD DUMMY 
//GO SYSIN DD ' 


= {TRK.(10.1)). 
= (CYL,(5,1)). 


C.7.1.3 NVT3DJCL3 


r 

u 
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C.7.2 Implementation Under VMBATCH 


NVT3D may be submitted under VMBATCH by submitting the file NVTBATCH. 


C .7.2.1 NVTBATCH EXEC 


&TRACE ALL 

CP LINK userid 191 100 RR password 
&IF &RETCODE NE 0 &GOTO -ERRORS 
inn r 

&IF &RETCODE NE 0 &GOTO -ERRORS 

GLOBAL TXTLIB VLNKMLIB VSF2FORT VPIUTIL CMSLIB TCSLIB AGIILIB PREVIEW 

GLOBAL LOADLIB VSF2LCAD 

&IF &RETCODE NE 0 &GOTO -ERRORS 

&IFN = datafile 

62 = DATA 

63 ^ A1 

64 = OUT 

FILEDEF 05 DISK &IFN &2 B 
FILEDEF 06 DISK &IFN &4 &3 

FILEDEF 08 DISK &IFN DISP &3 (RECFM U BLKSIZE 132 

FILEDEF 09 DISK &IFN STRAINX &3 (RECFM U BLKSIZE 132 

FILEDEF 12 DISK &IFN STRESS1 &3 (RECFM U BLKSIZE 132 

FILEDEF 11 DISK &IFN STRESSX &3 (RECFM U BLKSIZE 132 

FILEDEF 13 DISK &IFN ENERGY &3 (RECFM U BLKSIZE 132 

FILEDEF 14 DISK &IFN STRAIN 1 &3 (RECFM U BLKSIZE 132 

LOAD NVT3D (CLEAR START 

EXEC SENDFILE &IFN &4 &3 TO userid 

EXEC SENDFILE &IFN DISP &3 TO userid 

EXEC SENDFILE &IFN STRAINX &3 TO userid 

EXEC SENDFILE &IFN STRESS1 &3 TO userid 

EXEC SENDFILE &IFN STRESSX &3 TO userid 

EXEC SENDFILE SIFN ENERGY &3 TO userid 

EXEC SENDFILE &IFN STRAIN1 &3 TO userid 

&EXIT 

-ERRORS 

&TYPE EXEC TERMINATING BECAUSE OF ERRORS 

&TYPE RETCODE = &RETCODE 

&EXIT 

where userid is the account from which the program is being executed, password is the read 
password of the A-disk of the account, and datafile is the finite element data file. VMBATCH 
uses the compiled version of the source code; therefore, NVT3D TEXT must be on the A-disk 
in order for the program to run. Also, NVT3D must be compiled using Version 2 of VS Fortran. 
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The data file must be in the 80-column, fixed-format mode. When run under VMBATCH, NVT3D 
will not produce plots of deformed and undeformed meshes, and the job v.»ll abort if the plot 
flags are activated. 
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VIRGINIA TECH CENTER FOR 
COMPOSITE MATERIALS AND STRUCTURES 


The Center for Composite Materials and Structures 
is a coordinating organization for research and 
educational activity at Virginia Tech. The Center was 
formed in 1982 to encourage and promote continued 
advances in composite materials and composite 
structures. Those advances will be made from the 
base of individual accomplishments of the forty 
members who represent ten different departments 
in two colleges. 

The Center functions through an Administrative 
Board which is elected yearly and a Director who 
is elected for a three-year term. The general purposes 
of the Center include: 

• collection and dissemination of information 
about composites activities at Virginia Tech, 

• contact point for other organizations and 
individuals, 

• mechanism for collective educational and 
research pursuits, 

• forum and agency for internal interactions at 
Virginia Tech. 

The Center for Composite Materials and Structures 
is supported by a vigorous program of activity ;:t 
Virginia Tech that has developed since 1963. Research 
expenditures for investigation of composite materials 
and structures total well over seven million dollars 
with yearly expenditures presently approximating 


two million dollars. 

Research is conducted in a wide variety of areas 
including design and analysis of composite materials 
and composite structures, chemistry of materials and 
surfaces, characterization of material properties, 
development of new material systems, and relations 
between damage and response of composites. 
Extensive laboratories are available for mechanical 
testing, nondestructive testing and evaluation, stress 
analysis, polymer synthesis and characterization, 
material surface characterization, component 
fabrication, and other specialties. 

Educational activities include eight formal courses 
offeied at the undergraduate and graduate levels 
dealing with the physics, chemistry, mechanics, and 
design of composite fnaterials and structures. As of 
1984, some 43 Doctoral and 53 Master's students have 
completed graduate programs and several hundred 
Bachelor-ievel students have been trained in various 
aspects of composite materials and structures. A 
significant number cf graduates are now active in 
industry and government. 

Various Center faculty are internationally recog- 
nized for their leadership in composite materials and 
composite structures *hrough books, lectures, 
workshops, professional society activities, and 
research papers. 


Aerospace and Ocean 
Engineering 
Raphael T. Haftka 
Eric R. Johnson 
Rakesh L. Kapania 
Chemical Engineering 
Donald C. Baird 
Chemistry 
James E. McGrath 
Thomas C. Ward 
James P. Wightman 

Civil Engineering 
R. M. Barker 
Electrical Engineering 
loannis M. Besieris 
Richard O. Claus 


MEMBERS OF THE CENTER 


Engineering Science 
and Mechanics 
Hal F. Brinson 
Robert Czarr 
David Ditlarc 
Norman E. L ig 
John C Duke, ,.. 

Daniel Frederick 
O. Hayden Griftm, Jr. 
Zafer Gurdal 
Robert A. Heller 
Edmund G. Henneke, II 
Carl T. Herakovich 
Robert M. Jones 
liviu Librescu 
Alfred C. Loos 
Don H. Morris 
John Morton 
All H. Nayfeh 
Marek Pindera 
Daniel Post 


J. N. Reddy 
Kenneth L Reifsnider 

C. W. Smith 

Wayne W. Stine hcomb 
Surot Thangjitham 
Industrial Engineering and 
Operations Research 
Joel A. Nachlas 

Materials En jneering 

D. P. H. Hasselman 
Robert E. Swanson 

Mathematics 
Werner E. Kohler 
Mechanical Engineering 
Charles E. Knight 


Inquiries should be directed to: 

Center for Composite Materials and Structures 
College of Engineering 
Virginia Tech 
Blacksburg, VA 24061 
Phone: (703) 961-4969 
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